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INTRODUCTION 

A  goal  in  fan-out  filter  design  is  construction  of 
a  single  input,  denumerably  infinite  output  spectrum  analyzer 
that  exhibits  the  individual  terms  of  the  input  Fourier 
series.   This  entails  a  fan-out  filter  network  that  has  a 
resistive  input  driving  point  impedance  or  constant  inter- 
actance, and  each  fan-out  that  has  a  narrow  band-pass  filter 
transfer  function  and  predictable  linear  phase  shift. 

Though  the  spectrum  analyzer  is  the  final  goal,  this 
thesis  will  consider  a  simpler  problem  which  is  a  one  input- 
two  output  filter  designated  as  the  complementary  filter 
by  Zobel  (8).   Early  research  will  be  discussed  briefly 
with  special  attention  to  Fritzemeyer' s  (1)  conceptual  simpli- 
fications such  as  interactance  and  impedance  elision.   Con- 
stant-k  complementary  filter  designs  will  be  presented. 
Improvement  of  filter  interactance  and  transfer  function 
filter  characteristics  will  be  accomplished  by  the  aidentity 
algorithm  and  impedance  elision.   Fritzemeyer  investigated 
Newton  complementary  filters  and  these  are  not  considered 
in  this  thesis.  A  large  number  of  appendices  contain  proofs 
and  digital  computer  programs  necessary  for  the  development 
of  constant-k  complementary  filters. 

PREVIOUS  WORK 

General 

0.  J.  Zobel  (8)  was  concerned  with  improving  the  driving 


point  impedance  (DPI)  of  a  constant-k  filter  and  a  constant-k 
complementary  filter.   Zobel  added  series  and  shunt  elements 
in  front  of  the  constant-k  filter  and  then  empirically  deter- 
mined the  value  of  these  elements.   In  his  constant-k  com- 
plementary filter  the  shunt  annulling  branches  were  elimi- 
nated, similar  to  impedance  elision,  to  obtain  the  best  DPI 
results.   No  investigation  of  the  voltage  transfer  function 
(VTF)  was  made. 

Bode  (5)  was  primarily  interested  in  improving  the 
DPI  of  a  constant-k  filter  by  adding  a  conductance  control- 
ling network  and  a  susceptance  annulling  network.   Later 
Bode  (4-)  used  these  methods  on  constant-k  complementary 
filters  and  a  type  of  impedance  elision  with  pole-zero  ana- 
lysis to  improve  the  DPI. 

Guillemin  (10)  investigated  complementary  and  potent- 
ially complementary  constant-k  filters  with  added  series 
elements  that  had  properly  chosen  coefficients.   The  resulting 
filter  was  essentially  that  which  Zobel  and  Bode  had  sug- 
gested. 

Norton  (11)  working  with  complementary  filters  determined 
what  the  form  of  the  equivalent  DPI  needs  to  be  in  order 
to  exhibit  a  constant  interactance  and  then  attempted  to 
synthesize  the  complementary  filters  which  had  this  DPI. 
The  calculations  necessary  for  this  synthesis  are  rather 
complex  as  they  involve  the  solution  of  several  simultaneous 
nonlinear  equations. 


Szentirmai  (7)  suggested  DPI  improvement  by  adding 
resistive  pads  (networks)  in  each  fan-out  of  a  multiple 
output  filter.   This  had  the  disadvantage  of  high  losses 
due  to  the  pads,  although  it  does  make  each  individual  fan- 
out  insensitive  to  each  of  the  other  fan-out  networks. 

Rowlands  (6)  suggested  deletion  of  components  nearest 
the  paralleled  terminals  of  complementary  fan-in  filters. 
This  is  impedance  elision  alright  but  no  justification  is 
given. 

Fritzemeyer  (1)  used  these  early  efforts  and  King's 
(2)  approximate  identity  to  investigate  design  procedures 
for  obtaining  the  most  constant  DPI  in  constant-k  and  comple- 
mentary filters.   He  suggested  a  different  type  of  impedance 
elision  than  did  Bode  or  Zobel  and  his  suggested  aidentity 
algorithm  was  different  and  realizable  when  compared  with 
Norton's  method.   However  none  of  these  men  were  concerned 
with  the  VTF  of  their  resultant  filters  and  it  is  felt  that 
this  cannot  be  ignored.   Therefore  this  thesis  was  initiated 
and  completed  with  favorable  results  utilizing  the  suggestions 
of  the  referenced  authors. 

DEFINITIONS  AND  DESCRIPTIONS 

Lattice  Representation  of  Symmetric  Networks.   Any 

symmetric  network  may  be  represented  by  an  equivalent  lattice 

.  • ^- a  _, 0 

such  as:  ^-v^  ' — »    - 


or  the  couple,  -{a;b\,  where  a=f-,(s)  end  b=fp(s).   fa;b}will 
be  employed  in  place  of  the  lattice  in  textual  material 
for  compactness. 

Characteristic  Impedance.   Any  symmetric  network,  /a;b|, 
has  a  characteristic  impedance,  Z  ,  equal  to  /s/ab. 

Gonstant-K  Low-Pass  Filter.   A  symmetric  low-pass  filter 
can  have  a  characteristic  impedance  which  is  an  approximation 
of  a  desired  characteristic  impedance,  k  ohms.   For  instance 
the  prototype  low-pass  filter,  {Z-,/2;  (Z-j+^-Zp)/^,  has  a 
characteristic  impedance,  k^/l+CZ-j/^-Zp)  where  k=/s/z7zT, 
and  a  cut-off  frequency,  w  =2AyL7Cp~  if  Z,=L-,s  and  Zp=l/CpS. 
See  Fig.  1. 

Normalized  Low-Pass  Filter.   A  network  with  circuit 
components  adjusted  so  that  the  characteristic  impedance 
is  one  ohm  and  the  cut-off  frequency  is  one  radian  per  second 
and  is  designated  by  the  couple  (Z  ,  w  )=(1,1).   This  normal- 
ization  permits  computational  ease  without  lessening  general- 
ity.  It  is  always  possible  to  retrieve  specified  (Z  ,  w  ) 

o   c 

by  an  impedance  level  change  and  a  time  scale  change  (i.e., 
s  is  replaced  byTs). 

Constant-1  Low-Pass  Filter.   A  normalized  symmetric 
low-pass  filter  with  k=l  and  w  =1  rad./sec.   These  normali- 
zations  require  L-.=2  henrys  and  Cp*2  farads  for  the  prototype 
filter  of  Fig.  1  which  results  in  the  Constant-1  filter, 
/s;(s  +l)/s},  of  Fig.  -2.   Note  that  the  actual  characteristic 
impedance  is/Vl+s  and  this  is  an  approximation  to  k=l 

Zobel 's  M-Derivation.   Given  {a;b|  and  a  constant  m 


such  that  o<ra<l,  then  the  network,  /ma;  b/mj-  ,  is  a  network 
m-derived  from  the  given  network.   The  two  networks  have 
the  same  characteristic  impedance  and  the  order  of  cascading 
can  be  ignored.   Fig.  3  shows  m-derivation  with  both  a  T-net- 
work  and  its  lattice  equivalent. 

Network  Transposition.   The  transpose  of  a  four-terminal 
network  is  accomplished  by  transposing  the  network's  right 
and  left  terminal  pairs.   The  transposed  network  has  a  trans- 
fer matrix  with  the  a-,-,  and  app  elements  interchanged. 
If  network  transposition  leaves  the  network  unchanged,  then 
the  given  network  is  symmetric. 

Zobel  Filter.   A  symmetric  network  obtained  by  the 
process  of  m-deriving  a  T-network  prototype,  -fs;(s  +  l)/sV; 
bisecting  the  resulting  -TmsjCs  +  l)/msj';  and  inserting  the 
prototype  between  its  bisected  and  transposed  m-derivation. 
Fig.  4  illustrates  a  Zobel  filter  normalized  to  (1,1). 

Degenerate  Zobel  Filter.  A  symmetric  network  obtained 
by  taking  a  Zobel  filter  and  replacing  the  constant-1  filter 
with  the  identity  network.   Fig.  5  illustrates  a  degenerate 
Zobel  filter. 

Bisected  Zobel  Filter.   A  symmetric  network  obtained 
by  taking  a  Zobel  filter  and  replacing  the  constant-1  filter 
with  a  m-derived  7T~network«  Fig*  6  illustrates  a  bisected 
Zobel  filter. 

Aidentity.   A  ratio  of  polynomials  in  s  whose  numerator 
and  denominator  are  arranged  in  increasing  powers  of  s  and 
whose  ratio  approximates  the  constant,  1. 


Aidentity  Driving  Point  Impedance  (ADPI).   A  driving 
point  impedance  (DPI)  with  the  aidentity  property.   Such 
a  DPI  is  both  a  positive  real  function  (prf)  and  an  aident- 
ity. 

Aidentity  Order.   In  en  aidentity  if  the  first  p  suc- 
cessive  pairs  of  coefficients  of  s  ,  o<k<p-l,  in  the  numer- 
ator and  denominator  are  equal,  the  aidentity  order  is  p. 

Complementary  Filter.   A  low-pass  filter's  input  term- 
inals connected  in  parallel  with  a  high-pass  filter's  input 
terminals  results  in  a  complementary  filter.   The  cut-off 
frequencies  of  the  two  filters  are  the  seme. 

ASSERTIONS 

Assertion  1.   Two  different  symmetric  networks  can 
have  the  same  characteristic  impedance. 

Proof:   Given  |a;b>  and  /ma;b/m}.   Both  have  the  same 
Z   and  yet  they  have  different  circuit  components  and  transfer 
characteristics. 

Specific  examples  in  the  prototype  low-pass  filter 
context  are  given  in  Fi^s.  7,  8,  9,  end  10.   It  is  observed 
that  the  bisected  process  of  Fig.  10  yields  a  different 
example  network  having  the  same  Z   as  the  three  previous 
example  networks. 

Assertion  2.   Whereas  the  m-derived  networks  of  Assert- 
ion 1  have  e  Z   independent  of  ra,  the  ADPI's  of  the  m-derived 
networks  in  Assertions  3»  4,  and  5  are  dependent  on  m. 


Assertion  5.   The  Zobel  process  or  m-derived  half  end 
termination  process  improves  the  characteristic  impedance 
of  the  networks  of  Assertion  1  after  m  is  chosen  to  be  0.707. 
This  is  shown  in  Fig.  11. 

Assertion  4.   The  ADPI ' s  of  the  constant-1,  m-derived, 
composites  and  bisected  Zobel  filters  and  loads  are  different, 
but  each  has  the  same  aidentity  order.   In  contrast  to  Asser- 
tion 2,  the  ADPI ' s  orders  are  not  dependent  on  m  for  these 
cases  shown  in  Figs.  12,  13,  14,  15,  and  16. 

Assertion  *?.   The  aidentity  order  of  the  Zobel  filter 
and  the  degenerate  Zobel  filter  is  greater  by  two  than  the 
aidentity  order  of  the  Assertion  4-  filters;  and  in  contrast 
to  Assertion  4,  the  ADPI * s  orders  are  dependent  on  m=0.707 
for  the  cases  of  Fig.  17  and  Fig.  18. 

Conclusions 

The  Zobel  process  effects  can  be  summarized  as  follows: 

a.  When  the  characteristic  impedance,  Z  ,  is  dependent 
on  m,  there  is  a  resulting  improvement  of  Z   and 
ADPI  order  when  m  is  chosen  to  be  0.707. 

b.  When  the  characteristic  impedance,  Z  -  is  independent 
of  m,  there  is  no  improvement  of  Z   or  ADPI  order 
from  any  realizable  choice  of  m. 
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CONSTANT-K  COMPLEMENTARY  FILTERS 

Procedure 

In  the  previous  section  it  was  shown  that  the  Zobel 
process  improved  the  characteristic  impedance  and  ADPI  of 
low-pass  prototype  filters.   Now  we  shall  compare  the  fre- 
quency response  characteristics  of  the  constant-1  filter, 
Zobel  filter  (with  m  having  two  different  values),  constant-1 
complementary  filters  (modified  by  aidentity  algorithm) 
and  complementary  Zobel  filters  (modified  by  impedance  elision 
and  the  aidentity  algorithm).   The  filter  characteristics 
of  interest  are  interactance,  ADPI,  and  voltage  transfer 
function  (VTF).   Interactance  versus  frequency  is  a  unique 
filter  figure  of  merit  first  presented  by  Fritzemeyer  (1). 
Interactance  is  an  effective  filter  design  substitute  for 
characteristic  impedance.   The  interactance  freauency  res- 
ponse was  calculated  and  tabulated  on  the  IBM  1620  digital 
computer.   A  review  of  the  interactance,  the  computer  program, 
and  the  tabulated  input  and  output  data  are  given  in  Appendix 
C  and  Appendix  D.   The  ADPI  versus  frequency  was  also  computed 
on  the  1620  and  this  program  and  results  are  given  in  Appendix 
E.   Impedance  elision  is  presented  in  Appendix  A  and  the 
aidentity  algorithm  is  presented  in  Appendix  B.  All  of  these 
filters  have  the  ladder  configuration  and  a  computer  program 
was  written  to  determine  the  coefficients  of  each  filter's 
low-pass  VTF.   This  is  given  in  Appendix  F  with  the  computer 
program  and  the  resultant  data  given  in  Appendix  G.  Appendix 
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H  shows  the  computer  program  and  results  for  computing  the 
low-pass  VTF  frequency  response. 

Constant-1  Filter 

Considering  the  low-pass  constant-1  filter  of  Fig.  19 
first  will  give  a  basis  of  comparison  for  the  various  improve- 
ments of  the  filter  characteristics  of  the  Zobel  filters. 
The  ADPI  of  eauation  (1)  for  the  constant-1  filter  has  an 
aidentity  order  that  can  not  be  increased  by  the  aidentity 
algorithm  as  is  shown  in  Appendix  B.   Fig.  20  shows  the 
interactance  curve  which  increases  rapidly  between  0.3  rad/sec 
and  1.0  rad/sec. 

2    3 
A  +An  s+A0s  +Axs 

Zin  -   °  -   \     ?  <« 

B  +Bn  s+B0s 

O   1    <L 


where  A^  ■  1  B  ■  ■  1 

o  o 


kx   =  2  Bj  >  2        (2) 


A2  .  2  B2  •  2 

v  A   =2 
A^  -  d 

This  indicates  that  the  ADPI  aidentity  order  has  to  be 
increased  to  improve  the  interactance.   The  interactance 
is  verified  by  the  ADPI  curve.   The  low-pass  VTF  curve  does 
not  have  a  very  sharp  knee  and  the  VTF  phase  shift  curve 
is  linear  throughout  the  pass  band  as  shown  in  Fig.  21, 
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Zobel  Filter 

Using  the  common  value  of  0.6  for  m  in  the  Zobel  filter 
of  Fig,  22,  the  ADPI  of  equation  (3)  for  this  filter  has 
an  aidentity  order  of  one  and  has  the  characteristic  responses 
of  Fig.  23  and  Fig.  24. 


in 


2     3     4     5    6     7 
Co+ClS+C2   C3   G4S  *C5S  C6S  +C7S 

2    3    4-56 
I)o+I)ls+D2   D3    4s  +D5S  +D6S 


(3) 


where 


Co 

■ 

1 

cl 

= 

3. 

,84 

C2 

s 

5. 

.76 

C3 

= 

11, 

.63 

C4 

= 

6, 

.32 

C5 

s 

9. 

,1 

U6 

= 

1, 

.95 

7 


2.07 


Do  «  1 
Dx  =  3.2 


D. 


8.06 
6.64 


D4  «  10.81 
D^  -  3.07 
D6  =   3.26 


(4) 


However  when  the  aidentity  algorithm  is  applied  to  the  ADPI 
as  in  Assertion  5  on  page  7 »  equation  (5)»  which  has  an 
aidentity  order  of  five,  is  obtained  and  m  is  0.707. 


?    3    4    5    6    7 
Co+Cls+C2  +G3   C4S  +C5S  +C6S  +C7S 
Zin  =  £ 


(5) 


2     3     4     5    6 
DQ+D1s+D2s  +D^s^+D^s  +D,-S''+D£S 
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where 


Co 

= 

1 

cl 

s 

3.4-14 

C2 

= 

6.828 

C5 

s 

9.243 

C4 

■ 

8.992 

C5 

= 

6.682 

C6 

■ 

2.914 

C7 

= 

1.457 

Do  - 

1 

Dl  - 

3.414 

D^  ■ 

6.828 

D3- 

9.243 

D4  " 

8.992 

D5» 

5.828 

D6  - 

2.914 

(6) 


This  results  in  a  Zobel  filter  with  constants  as  shown  in 
Fig.  25  and  characteristic  responses  of  Fig.  26  and  Fig.  27. 

In  comparing  the  characteristic  responses  for  the  two 
choices  of  m,  it  is  obvious  that  the  interactance  for  m=0.707 
is  much  better  during  the  pass-band  than  for  m=0.6.   This 
was  expected  because  of  the  difference  in  the  aidentity 
order  of  the  two  ADPI'S.  It  is  difficult  to  make  a  conclusion 
about  the  interactance  improvement  of  the  Zobel  filter  over 
the  constent-1  filter  because  of  the  extreme  variations 
near  w  for  the  Zobel  filter. 

It  is  worth  noting  that  Fritzemeyer's  (1)  interactance 
curve  for  the  m-derived  filter,  with  m=0.6  and  ADPI  aidentity 
order  of  3»  is  essentially  the  same  as  the  interactance  curve 
of  Fig.  26  for  a  Zobel  filter  with  m=0.707  and  ADPI  aidentity 
order  of  5.   The  validity  of  Fritzemeyer's  curve  is  auestioned. 

It  should  be  noted  that  the  low-pass  VTF  response  for 
m=0.6  has  a  slightly  sharper  knee  than  the  one  for  m=0.707. 
However  both  have  much  better  low-pass  VTF  responses  than 
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the  constant-1  filter  and  both  are  approximately  0.35  at 

\v  .   The  VTF  phase  shift  curves  for  both  values  of  m  are 
c 

not  as  linear  as  was  the  constant-1  phase  shift  curve.   The 
phase  shift  is  rather  linear  from  zero  to  about  0.75  rad/sec, 
then  the  curve  becomes  slightly  nonlinear.   Also  the  Zobel 
filter  phase  shift  slope  is  steeper  than  the  slope  for  the 
constant-1  filter.   The  ADPI  phase  shift  is  improved  consid- 
erably over  the'm»0.6  Zobel  filter  and  the  constant-1  filter. 

Complementary  Filters 

Having  looked  at  the  characteristics  of  the  prototype 
low-pass  constant-1  and  Zobel  filters,  we  will  now  observe 
the  characteristics  of  complementary  filters  when  the  aident- 
ity  algorithm  and  impedance  elision  are  applied. 

Type  1  Constant-1  Complementary  Filter 

When  two  constant-1  complementary  filters  are  connected 
as  shown  in  Fig.  28,  the  ADPI  of  eauation  (7)  results. 
No  aidentity  algorithm  is  possible  with  this  ADPI  as  all 
components  have  fixed  values.   The  symmetry  of  numerator 
and  denominator  coefficients  should  be  observed.   This  is 
expected  because  of  the  mathematical  operations  in  determining 
the  eauivalent  parallel  impedance  when  high-pass  and  low-pass 
impedance  functions  are  involved.   Appendix  I  provides  a 
further  explanation  of  this  coefficient  symmetry.   Inspection 
of  the  interactance  curve  of  Fig.  29  reveals  that  the  comple- 
mentary nature  of  this  filter  will  result  in  a  symmetrical 
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curve   about   w    . 
c 


'cf 


E  +E1s+E2s   +E5s^+EZfs   +EcS^+E6s 
F0+F1S+F2S  -»-F,s^+F^s4+Fc8^+F6s 


(7) 


where 


Eo 

s 

E6 

s 

2 

El 

s 

E5 

s 

6 

E2. 

= 

E4 

= 

10 

E3 

■ 

13 

F- 


r4 
18 


2 

8 
16 


(8) 


Note  that  although  the  interactance  curve  is  not  flat,  it 

does  not  have  the  extreme  variations  about  w„  that  the  con- 

c 

stant-1  filter  and  Zobel  filter  each  have.   Its  VTF  response 
of  the  low-pass  fan-out  is  better  than  the  constant-1  filter, 
but  not  quite  as  good  as  the  Zobel  filter.   It  is  approx- 
imately 0.3  at  w  and  doesn't  reach  0.1  until  w=1.5  rad/sec. 
However  its  low-pass  VTF  phase  shift  curve  is  as  linear, 
but  with  less  slope,  as  the  VTF  phase  shift  of  the  Zobel 
filter.   Both  have  the  same  type  of  nonlinearity  at  w  . 
The  ADPI  phase  shift  response  has  deteriorated  by  adding 
the  complementary  high-pass  filter  in  parallel  with  the 
constant-1  filter. 


Type  2  Constant-1  Complementary  Filter  with 

Parameters  C  and  C ' 


To  improve  the  interactance  of  the  complementary  con- 
figuration, parameters  C  and  C',  were  added  to  the  reactances 
nearest  the  input  terminals  of  each  filter  (Fig.  3D  so 
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that  the  aidentity  algorithm  could  be  used  on  the  ADPI  of 
eouation  (9).   Again  as  stated  in  the  previous  complementary 
filter,  G  =G6,  G^Gj-,  G2=G4,  H0sH6i  H^Hr,  and  H2«H4.   The 
aidentity  algorithm  is  next  applied  which  set  G-,»H,  and 


Gc=Hr. 


P     3    4     S    6 
GA+Gn  s+G0s-+Gzs^+G,,s  +Gcs''+Gc.s 

Z„r   =  — S— i § 5 2 § § (9) 


'cf 


Ho+giS4-H2s2+H$s5+H4s4+H5s5+H6s6 


where       G^  =  2C •  =  Gc  =  2C 

O  D 

G1  =  2C'(2  +  C)  =  G^  =  2C(2  +  C) 

G2  «  1  +  3C  +  6CC  =  G^  =  1  +  3C  +  6CC 

G*  -  2  +  C  +  C  +  9CC  (10) 

H  «  2C'  =  Hc  »  2C 

O  D 

Hx  =  2(1  +  3C)  =  H^  =  2(1  +  3C) 

H2  -  5  +  2C  +  9C  =  H4  =  5  +  2C(  +  9C 

H^  =  6(1  +  C  +  C) 

These  two  simultaneous  equations  yield  the  results,  C«C'«1.618 
and  a  second  order  ADPI.   This  suggests  that  the  dual  com- 
ponents of  the  complementary  filters  should  h^ve  the  same 
coefficients.   This  observation  will  be  used  on  later  com- 
plementary configurations.   Using  this  value  for  C  and  C, 
Z  -  has  the  coefficients  of  equation  (11)  and  the  frequency 
response  characteristics  for  the  complementary  filter  of 
Fig.  31  given  in  Fig.  32  and  Fig.  33.   By  adding  the  para- 
meters C  and  C  and  using  the  aidentity  algorithm  on  Z  -, 
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GQ  -  G6  -   5.24 
Gx  .  Q     m   11.7 


0 


2  -  G4  -  21.55 


Ho 

■ 

H6 

m 

5. 

,24 

Hl 

" 

H5 

a 

11, 

.7 

H2 

- 

H, 

s 

22, 

.8 

H3 

2 

25. 

,42 

G^  =  28.79 

the  interactance  of  Fig.  32  results  which  during  the  larger 
part  of  the  high  and  low  pass-bands  is  much  better  than  for 
Type  1.   However  the  interactance  response  about  w  has 
deteriorated  from  Type  1.   The  VTF  response  of  the  low-pass 
fan-out  has  been  improved  from  Type  1  low-pass  as  it  is 
flat  for  more  of  the  pass-bend,  then  has  a  sharper  knee, 
a  0.312  value  at  w  ,  and  a  0.1  value  at  w=1.25  rad/sec. 
The  low-pass  VTF  response  is  similar  to  that  of  the  Zobel 
filter  with  m=0.707.   The  low-pass  VTF  phase  shift  curve 
is  nearly  identical  with  that  of  Type  1,  with  the  same  non- 
linearity  at  w  .   The  ADPI  phase  shift  response  during  most 
of  the  pass-bands  is  much  better  than  for  Type  1  except  for 
the  larger  variations  about  w  .   It  does  have  the  same  gen- 
eral  shape  as  the  ADPI  phase  shift  curve  of  the  Zobel  filter. 

Type  3  Oonstant-1  Complementary  Filter  with 

Parameters  C  and  C, 

Using  the  idea  of  Type  2  that  dual  components  of  the 
complementary  filter  should  have  equal  coefficients,  the 
parameters  C  and  G-,  were  added  to  the  constant-1  complementary 
filter  (Fig.  34)  with  the  idea  of  using  x;he  aidentity  algorithm 
to  improve  the  interactance  and  ADPI.   The  ADPI  of  equation  (12) 
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I  +I,S+I~S  +  I2S  +I..S  +IC8^+ICS 

Z  -   -  — £— i § 2 2 1 § (12) 

?    3    4    5    6 

J  +J-.  S+JpS  "+J,8y+J^S  +JrS  +«L-s 


where  I  -  Ig  -  CC1  «  3.732 

Ix  .  I5  -  0(^(2  +  c)  =  13.93 

12  -  I4  -  1  +  C  +  CC1(1  +  C  +  CC^)   =  26.86 

13  =  2(1  +  C)  +  C2(l  +  2C12)  «  36.33  (13) 

Jo  =  J6  *  CC1  =   3*752 

Jx  .  j_  ■  c1(l  +  C  +  CC1)  -  13.93 

J2  =  j^  =  1  +  c  +  2C1  +  CC1(1  +  2C1)«  26.86 

J5  =  2[l  +  C1  +  CC1(1  +  C^)]        =  29.86 

results  from  the  circuit  of  Fig.  34.   Using  the  aidentity 

algorithm  on  the  simultaneous  equations,  I  «J^,  I-,  -Jt  ,  and 
"  T   o  o'   1  1 ' 

Ip=Jo,  yields  C=1.732  and  0^=2.155.  Comparing  the  resultant 
component  coefficients  of  Type  3  with  Type  2,  indicates 
that  the  changes  are  rather  small.   This  implies  that  the 
interactance  may  not  be  improved  much  by  this  approach, 
although  the  ADPI  order  is  now  three.   Examination  of  this 
filter's  frequency  response  characteristics  given  in  Fig.  35 
and  Fig.  36  confirms  the  above  suspicions.   There  is  even 
a  more  extreme  interactance  variation  about  w  for  Type  3 
than  for  Type  2  which  agrees  with  the  large  difference  between 
I,  and  J-,.   The  low-pass  VTF  response  even  has  a  more  rounded 
knee  than  did  Type  2,  however  it  reaches  0.1  at  w=1.15  rad/sec 
instead  of  w=1.25  rad/sec  as  did  Type  2.   The  one  bright 
spot  is  the  improvement  of  linearity  in  the  low-pass  VTF 
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phase  shift  curve  about  w  which  is  not  present  in  Type  2. 
Along  with  the  other  deteriorating  features,  the  ADPI  phase 
shift  curve  has  greater  variations  about  w  than  did  Type  2, 
although  it  has  less  phase  shift  during  most  of  its  pass-bend 
than  Type  2.   It  seems  that  by  increasing  the  aidentity 
order  by  one  there  is  only  a  small  improvement  in  the  phase 
shift  response. 

Type  4  Constant-1  Complementary  Filter  with 
Parameters  C,  C-,  and  Cp 

The  results  of  Type  3  seem  to  imply  that  any  further 
attempt  at  improving  the  ADPI  and  interactance  by  addition 
of  a  third  parameter  would  be  futile.   However  when  the 
three  component  parameters  were  used  as  shown  in  Fig.  37* 
the  ADPI  of  eauation  (14)  resulted.   The  aidentity  algorithm 
using  K  =L  ,  K.,=L..  ,  Kp=Lp,  and  K;z=L;r  yields  C=1.5»  0,-1. 33» 
and  Cp=0.5  which  then  results  in  an  aidentity  order  of  seven. 
This  bonus  in  aidentity  order,  of  course,  comes  from  the 
property  of  coefficient  symmetry  of  the  ADPI.   The  frequency 
response  characteristics  of  Fig.  38  and  Fig.  39  reflect 
an  ADPI  of  one  for  all  freauencies  and  a  constant  interactance, 

K  +K,  s+K^+K^+K.s^+Kc-s^+K^s6 
Z    .   =        °      X § 2 H 2 § (14) 

2  5  4  5  6 

L  +Ln  s+L0s  +L-,s^+L.. s  +Lcs^+Iv-s 
o      1        2  3  4  ^  6 


27 


where 


1C      =   K,-    «  CC-,Co 
o  o  Id 


s 

0.998 

■ 

3.99 

= 

7.98 

■ 

9.95 

= 

0.998 

= 

3.99 

= 

7.98 

■ 

9.95 

K1    =    K^   =   CC^l   +   CC2   +   C^2) 

Kp   -  K^   =  C   +   C2  +   0(^(0   +   C2  +  CC^Cg) 

K,   =    1    +   C2  +   C22  +   2CC2  +   C2C12(1   +  C22) 

o  b  Id 

i^  =  L.  =  o1c2  +  cc^i  +  c1c2) 

L2  =   L^   =  C    +   C-f    +   C2  +  C1C2      +  CC1(C1 

+  c2  +  C^Cp) 
L,   =   2[l   +  C^g  +  0(^(1   +  C^Og)] 

This  is  the  ultimate  in  driving  point  impedance  character- 
istics for  most  filters.   The  low-pass  VTF  response  has 
been  improved  by  the  w  value  being  0.353  while  the  rest 
of  the  pass-band  response  remains  very  nearly  the  same  as 
Type  3-  This  finding  is  in  disagreement  with  what  Fritze- 
meyer  (1)  stated  would  be  the  response.   The  ADPI  phase 
shift  is  a  straight  line  at  zero  degrees  as  is  expected. 
The  low-pass  VTF  phase  shift  has  less  nonlinearity  about 
w  than  did  Type  3  and  its  slope  is  also  less,  which  gives 
it  a  very  good  linear  phase  shift  response. 

Type  5  Complementary  Zobel  Filter 

As  has  been  the  procedure  of  the  two  previous  sections, 
the  Zobel  filter  configuration  is  used  to  improve  the  ADPI 
of  the  complementary  filter.   Using  the  network  of  Fig.  40 
results  in  the  ADPI  of  equation  (16).   The  determination 
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Z„,  -  -2=2 (16) 


'of 
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Ev" 


n=o 


where 

a  =   1   +  m 

b  =   1   -  in 


4^2 


Mo   =   M14   =  2a-b 


Ml  =  M13  =  ?a3b(1+2a2b) 

M2  =   M12   =   2a5b  [2(a-b)   +   4a2b] 

M     =  Mn   =   4aS[a2b(l-2b)   +   4  +   b]    +    a2[2  +  b(2+a)] 

M4  =  M1Q  =   2a2{a[2  +   b(2+a)]  (l+a3b5)+(l+2b)(l+2a5b) 

+   4a2b(2+b)} 
Mc  =   M9  »   2a5{[2  +   b(2+a)](2+a2b)   +    2(l+2b)[l  (1?) 

+   a2b(2+b)]    +   2a4b5(2+b)}  +   4a 
M6  =  Mg  »   2a4[2  +   b(2+a)][l   +   2b  +   ab(2+b)]    +   4a5[2 

+  4b  +  a5b2(2+b)  +  a4b5]  +  8a2  +  2a 


M 


7 


=  1  +  20a2  +  4a6b2[a2b2  +  1  +  (2+b)2  ]  +  a4{[2 
+  b(2+a)]2  +  4(l+2b)2} 


N  *  N,..  =  2a4b2 

O      14 


Nx   =   N15   =  4a5b(l+a2b) 

N2  =  N12  =   2a2{1  +   ab[2  +  b  +   4a(l+ab)]} 
N3  =  Nll   =  2&2b  +  b(2+a)]  +  4a5{i  +   ab[6  +  b 
ib(l+a)]} 


+   ai 
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N4  *  N10  *  (1+2a3b)[2a(1+a)  +  2a2(l+2b)]  +  8a3[l 

+  ab(2+b)]  +  [2  +  b(2+a)][2a5(2+a5b2)] 
H     =  N9  -  8a  +  4a2(l+a)[l  +  a2b(2+b)]  +  4a5{(l+2b)(l+a) 

+  [2  +  b(2+a)]  (2+a2b)  +  a5b2} 
N6  =  Ng  =  4a  +  16a2  +  2a5  [2  +  b(2+a)j  [l  +  a  +  a2b(2+b)  (17) 
+  2a(l+b)]  +  8a2[l  +   +  a(l+2b)]  +  4a5b[l 
+  s2b2  +  ab(2+b)] 
N7  =  2  +  40a2  +-8a5((l+a)(l+2b)  ♦  a2b(ab  +  2  +  b)] 
+  2a4  [2  +  b(2+a)] 

of  the  ADPI  coefficients  in  algebraic  form  was  quite  difficult. 
It  required  the  product  of  two  seventh  order  polynomials 
with  algebraic  coefficients.   Because  of  this  tedious  task 
a  computer  program  wan  written  to  eliminate  the  cataloging 
and  tabulating  errors  that  developed  when  the  multiplying 
was  done  by  hand.  The  program  is  illustrated  in  Appendix  J. 
When  the  aidentity  algorithm  is  applied  to  this  ADPI,  M  =N 
gives  no  information  as  they  are  identical  and  M,=N,  implies 
that  0=1  for  all  values  of  m  so  that  equation  (16)  will 
be  a  first  order  ADPI  for  all  realizable  values  of  m.   Fritze- 
meyer  (1)  assumed  m=0.6  and  obtained  the  interactance  response 
of  Fig.  41.   From  the  previous  sections  it  seems  that  the 
m=0.707  would  produce  a  better  interactance  response,  but 
the  ADPI  being  limited  to  a  first  order  aidentity  gives 
very  little  promise  of  much  improvement.   The  aidentity 
algorithm  could  be  used  at  this  point,  but  Fritzemeyer  (1) 
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suggested  an  additional  process  for  complementary  filter 
characteristic  response  improvement  which  will  be  called 
impedance  elision  as  detailed  in  Appendix  A.   Impedance  eli> 
sion  will  yield  the  degenerate  complementary  configuration. 


Type  6  Degenerate  Complementary  Zobel  Filter  with 
Parameters  C  and  m 


Applying  impedance  elision  to  the  Type  5  Filter  and 
adding  the  parameter  C  with  m  results  in  the  complementary 
configuration  of  Fig,  4-2.   The  two  parameters  allow  the  aident- 
ity  algorithm  to  be  used  on  the  ADPI  of  eauation  (18). 

All  values  of  C>  eaual  P^  because  both  are  identical  express- 

oo  e 

ions.      When  0,    is   equated  with  P,  ,    (l+m)C   =   (l+/s/5)/2  results 

p 

and  when  0~  is  equated  with  P~,  (l+m)C  =  l/(l-m  )  results. 

These  two  simultaneous  equations  result  in  m=0. 618033988 
(Golden  Mean)  and  C=l  and  a  third  order  ADPI.   The  numerical 
coefficients  of  equations  (20)  result  when  the  determined 
values  of  m  and  C  are  substituted  into  the  algebraic  coeffic- 
ients of  equations  (19).   It  should  be  noted  that  impedance 
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Z  vn 


Zcf  =  -S=2 (is) 

CI     10 


Zvn 


n=o 
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where 

a   ■   1   +   m 

b   =   1   -  m 

0     -  0,~   ■   2Ca*b 
o         io 


01   =   09   =   2Ca2[l  +   a2b(l+C)] 


°2 


=   0g  =   a2{b  +  C(2+b)    +   2Ca[l  +  C   +    ab(l+C+Ca)]J 
0,   =   0n  =   a(l+2Ca2)(l+C+Ca)   +   2Ca5b(b+Cb+2C)   +   a3(l 

+C)(b-f2C+Cb) 
04  =  06   =  4C2a5b  +   2Ca4(b+Cb+2C)   +   a2(l+C+Ca)  [ab 

+   aC(2+b)   +    1  +  CJ  +   a  +   aC 
0     *  4C2a4(l+a2b2)   +   a4(b+bC+2C)(b+2C+Cb)   +   a2(l+C 
+Ca)2  +    a2(l+C)2  +    1 

Po   =   P10   =   2Ca5b  (19) 

P1   =  P9  =   2a2{b  +  c[l  +   ab(l+a)]} 

P2  =  Pg  =  2a  [l  +   a2b(l+C)]    +   2Ca2[l  +   a  +   a2b(2+b)] 

+   a2(b+2C+Cb) 
P_   xx  p_   »    (l+C+Ca)(a+2a5b)   +   2a2(l+C)   +   2a  +   ab 

+   2Ca5(2+b+2ab)   +   a2(l+a)(b+2C+Cb) 
P^  s  P,    s   2a4b(b+Cb+2C)    +    (l+C+Ca)(a+3a2)   +   a2(2 

+b)[l  +  C    +    a(b+2G+Cb)]  +   4Ca5(l+a2b2) 

+   2a  +  Ca  +    1 
P5  =  2[l  +   a(l+C)(l+e)   +   a2(2+b)(l+C+Ca)   +  4Ca4b 
+   2a3(b+Cb+20)] 
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Po  ■  P10 

" 

3.236 

Pl  "  P4 

s 

15.708 

P2  "  P8 

■ 

43.125  (20) 

P?  =  P7 

= 

81.485 

P*  ■  p6 

B 

118.374 

P5  = 

131.846 

0Q  »  010  -  3.236 

01  =  0^  =  15.708 

02  =  08  =  43.125 
0$  =  0?  =  83.339 

04  =  06  =  119.374 

05  =  135.992 

elision  and  added  C  parameter  have  increased  the  aidentity 
order  by  two.   The  frequency  response  characteristics  ere 
given  in  Fig.  43  and  Fig.  44  for  this  degenerate  complemen- 
tary configuration.   The  interactance  has  been  greatly  im- 
proved by  impedance  elision  and  the  aidentity  algorithm  as 
compared  with  Type  5»   This  filter  has  by  far  the  best  low- 
pass  VTF  response  of  any  type  presented.  At  w  it  has  a 
value  of  0.32  and  reaches  0.1  at  w=1.08  rad/sec.   The  low-pass 
VTF  phase  shift  is  nearly  linear  to  approximately  0.8  rad/sec 

but  then  displays  the  usual  nonlinearities  about  w  .   The 

c 

ADPI  phase  shift  is  nearly  constant  from  zero  degrees  until 
0.8  rad/sec  but  then  becomes  nonlinear  about  w  as  has  been 
the  character  in  the  previous  types.   These  low-pass  char- 
acteristics substantiate  the  thinking  behind  the  use  of 
impedance  elision  as  they  are  very  nearly  the  same  as  those 
of  Fig.  16  and  Fig.  17  for  the  Zobel  filter  of  Fig.  15. 
Also  the  high-pass  series  arm  impedance  nearest  the  input 
terminals  of  Type  6  is  1.618/s  where  the  shunt  capacitive 
impedance  of  Fig.  15  is  1.414/s.  And  the  high-pass  shunt 
arm  impedance  nearest  the  input  terminals  of  Type  6  is  s/2 
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where  the  shunt  inductive  impedance  of  Fig.  15  is  s/1.414. 
This  implies  that  the  first  series  impedance  and  shunt  imped- 
ance of  the  high-pass  complementary  filter  of  the  degenerate 
Zobel  filter  takes  the  place  of  the  deleted  branch  in  the 
low-pass  complementary  filter.   The  characteristic  responses 
of  this  elision  or  sluring  have  been  shown  to  oe  quite  favor- 
able and  that  it  is  these  first  two  impedances  that  appear 
to  have  the  most  effect  on  the  other  complementary  filter. 

Type  7  Degenerate  Complementary  Zobel  Filter  with 
Parameters  C,  C-,  and  m 

An  additional  attempt  was  made  to  improve  the  ADPI 

of  the  Type  6  degenerate  complementary  configuration  by  the 

addition  of  another  circuit  parameter,  C-,,  as  is  shown  in 

Fig.  4-5.   From  this  network  the  ADPI  of  equation  (21)  is 

obtained  with  the  algebraic  coefficients  of  equations  (22). 

Obtaining  these  coefficients  required  a  considerable  amount 

of  toil  similar  to  that  experienced  with  Type  6.   Applying 

the  aidentity  algorithm  to  these  coefficients  results  in 

Q  =  R  which  gives  no  information  as  both  expressions  are 
^o    0 

identical.   Since  there  are  three  parameters,  at  least  three 
other  sets  of  eauations  (22)  coefficients  will  have  to  be 
equated.   The  result  of  Q,  eauated  with  R,  is  equation  (23) » 
of  Qp  equated  with  R~  is  eauation  (24),  and  of  Q,  eauated 
with  R,  is  equation  (25). 
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where 

a  =   1   +   m 

b  =  1  -  m 


3i 


Qo   -   Q10   =   a-bCC1 

Q1   »  Q     =   a^bCOj^Cl+C)   +   a2^ 

Q2  =  Q8  =   a^bc[c   +  C1(l+Cm)]    +   a5CC1(l+C)   +   a^jc 

+   bC1(l+C)] 
Q     ■  Q„  «   a5bc[c   +   bC1(l+C)]    +   a5c[c   +  C1(l+Cm)] 

+   a^CLfC)^   +   b01(l+C)]   +   aC^C   +  0^(1 

+  Cm)]  (22) 

Q4  =  Q6   =   a5bC12  +   a4c[b   +   bC1(l+C)]    +   a5[c   +   bC1(l 

+C)][c   +  C^l+Cm)]    +   a2C1(l+C)[c   +  0^(1 

+Cm)]   +   aC12(l+C) 
Qc  =   a6b2C2  +   a4C2  +   a4[c   +   bC]L(l+C)]2  +   a2[c   +  0^(1 

+Cm)]2  +   a2C12(l+C)2  +  Cj2 
Ro   =  H10   *   a5bCC1 

R.    ■  Rq  ■  a;?bC(l+C1m)  +   a  CC1  +   a  bC1 
R2  =  R8  =  a^bC(l+bC1)   +   a^Cd+Ojm)   +   a2C1[c  +  bC1(l 

+C)]   +   a5bC1(l+C)   +   aC1 
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R3  =  R7 


-   a4bC   +    a5C(l+bC1)   +    a^l+C^m)  [c   +   bC^d+C)] 
+    aC-^C   +  C1(l+Cm)]+   a5b[c   +   C1(l+Cm)] 
+    a2C1(l+C)   +   aC1(l+bC1) 
R^  =  R6   =   a5b2C   +   a5C   +    a5(l+bC1)[c   +   bC^l+C)]  (22) 

+   aCl+C^)^   +  C1(l+Cm)]    +    aC12(l+C) 
+    a4b[c   +   bG1(l-fC)]    +    a2[c   +  C]L(1+Cm)] 
+   a2C1(l+C)(l+bC1)   +  C1(l+C1m) 
R5  -  2{a4bC   +    a5[c   +   bC1(l+C)]   +    a2(l+bC1)[c   +  (^(1 
+Cm)]    +   aC1(l+C)(l+C1m)   +  C^2  !■ 


aC(C1+aCC1-l)    -  C1    =  0  (23) 

a5bC  [c-l+C1m(l+G)]+a2C1(l+G)(C-b)-.aC(l+C1m)-C1   «  0      (24) 

a4bc[p   +  bC1(l+C)]+   a5C  [cCj   -  bd+C^]    +   a2{cc1[2 

+   1*^(2+0  )]    +   b(C12-C1-C)    -  C   +  C2J+   a{c  (25) 

+  C1[b(l+C1m+CC1m)    -   l]  }-  C1(l+bC1)    *  0 

The  method  of  solving  these  three  simultaneous  equations 
is  rather  laborious  so  it  has  been  placed  in  Appendix  K. 
Their  solutions  give  the  following  values  for  C,  C-,  ,  and  m: 

C  =»  0.92041 
Cx  =  0.51437 
m  =  0.71433 

With  these  quantities  and  the  computer  program  of  Appendix  C, 
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the  ADPI  coefficients  of  equations  (22)   were  evaluated  and 
these  results  are  given  in  equations  (26).   This  shows  that 


%  •  Qio 

=   0.6814 

Ro  a  R10 

x   0.6814 

<?i  -  Q9 

-  3.6347 

Rl  =  R9 

-   3.6347 

Q2  =  Q8 

=  10.4251 

R2  a  R8 

-  10.4251 

(26) 

Q5  ■  Q7 

=  20.4523 

R^  ■  Ry 

=  20.4523 

Q4  =  Q6 

=  29.9044 

H4  -  R6 

»  29.9370 

Q5  = 

33.9341 

R5. 

33.7731 

the  ADPI  aidentity  order  has  been  increased  by  one  over 
Type  6  to  a  fourth  order.  Examination  of  Q^,  and  R^  reveals 
that  it  is  very  nearly  a  fifth  order  aidentity.  The  filter 
frequency  response  characteristics  are  given  in  Pig.  46 
and  Pig.  47.   The  small  improvement  of  the  frequency  response 
characteristics  of  this  type  over  Type  6  is  very  similar 
to  the  small  improvement  that  Type  3  had  over  Type  2.   In 
both  of  these  instances  the  added  parameter  was  in  the  shunt 
arms.   This  would  seem  to  indicate  that  the  shunt  arms  do 
not  have  as  much  influence  on  the  frequency  response  char- 
acteristics as  do  the  series  arms.  As  in  Type  4  when  the 
additional  parameter  was  placed  in  the  remaining  series 
arm  the  characteristics  improved  quite  significantly,  it 
is  believed  that  the  same  improvements  would  be  observed 
in  this  filter.   However  it  should  be  noted  that  solving 
for  the  four  parameters  would  be  a  tremendous  task  beyond 
the  scope  of  this  thesis.   The  frequency  response  charact- 
eristics of  this  filter  contradict  Fritzemeyer's  (1) 
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Fig.    19.      Constsnt-1   filter. 
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Fig.    22.      Zobel   filter  with  m=0.6. 
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Fig.  23.   Zobel  filter  characteristics  with  m=0.6. 
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Fig.  25.   Zobel  filter  with  m=0.707. 
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Fig.  28.   Type  1  constant-1  complementary  filter. 
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Fig.  29.   Type  1  freouency  response  characteristics. 


40' 


0' 


H 

-40" 

A 

S 

E 

-80° 

S 

-120° 

H 

I 

-160° 

F 

T 

-200° 

-220' 


11—  ARG  Z(jw) 

o-ARG  TQw) 


J L 


I      » 


0.0 


1.0 


0.2     0.4    0.6    0.8 
w(radians  per  second) 
Fig.  30.   Type  1  frequency  response  characteristics. 


1.2 


43 


in. 


l-O- 


Fig.  31.   Type  2  constant-1  complementary  filter  with  para- 
meter C  and  C '  . 
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Fig.  32.   Type  2  frequency  response  characteristics. 
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Fig.  34.   Type  3  constant-1  complementary  filter  with  para- 
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Fig.  37.   Type  4  constant-1  complementary  filter  with  para- 
meters C,  C-,  and  Cp. 
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38.   Type  4  freauency  response  characteristics. 
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Fig.  39.   Type  4  frequency  response  characteristics. 
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Fig.  40.   Type  5  Complimentary  Zobel  filter. 


48 


A 


10. 

6. 

4. 

2. 

1.0 

0.6 
0.4 

0.2 
0.1 


1 


. I 


t  i  i 


■  ■  I  ...  I 


0.1   0.2  0.4  0.6   1.0   2.     4.  6.   10 

w  (radians  per  second) 

Fig.  41.   Interactance  of  Type  5  complementary  Zobel  filter 
with  m=0.6. 

C(m+l)s    s ms 


-T-1/2S 


(1-ni  )s/ra 


1/ms 


in. 


C(m+l)/s    1/s 


m/s 


ib 


If— -^df-x 


s/2 


(1-ni  )/ms 


s/m 


in. 


Fig.  42.   Type  6  degenerate  complementary  Zobel  filter  with 
parameters  C  and  m. 
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Fig..  45-   Type  7  degenerate  complementary  Zobel  filter  with 
parameters  C,  C-,  and  m. 
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Fig.  47.   Type  7  frequency  response  characteristics. 
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assertions  as  they  have  not  deteriorated  from  those  of  Type  6, 
but  have  actually  been  improved. 

CONCLUSIONS 

Using  the  procedures  of  Zobel,  Bode,  Guillemin,  Norton, 
Rowlands,  and  Fritzeraeyer  the  Zobel  process,  aidentity  algor- 
ithm, and  impedance  elision  are  defined,  codified,  end  veri- 
fied as  methods 'of  improving*  the  frequency  response  charact- 
eristics of  constant-k  and  constant-k  complementary  filters. 
It  is  shown  that  these  methods  not  only  improve  the  ADPI 
or  interactance,  but  also  improve  the  VTF  and  linearize 
the  phase  shift. 

In  the  course  of  this  verification  several  related 
topics  of  interest  were  studied  and  are  included  in  the 
appendices.   These  topics  are  the  FORGO  program  for  eval- 
uating the  coefficients  of  a  pseudo-generalized  "ladder" 
VTF,  the  FORTRAN  program  for  coefficient  evaluation  of  the 
product  of  two  polynomials  with  algebraic  coefficients, 
the  proof  of  the  complementary  filter  ADPI  coefficient  sym- 
metry, and  the  ADPI  parameter  evaluation  for  a  Type  7  filter. 
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APPENDIX  A 

Impedance  Elision.   Webster's  definition  of  elision: 
A  cutting  off,  especially  of  a  vowel,  for  the  sake  of  meter 
or  euphony;  the  dropping  or  partial  pronounciation  of  a 
final  vowel  before  an  initial  vowel  in  the  next  word.   Fritz- 
emeyer  (1)  used  impedance  elision  without  naming  it.   Given 
the  complementary  filter  configuration  of  Fig.  48,  it  was 
noted  thst  if  the  two  shunt  arms  nearest  the  input  terminals 
were  removed  that  the  effects  of  each  could  be  substituted 
by  the  first  series  component  and  the  next  shunt  arm  element. 
Or  when  the  shunt  arms,  ^rjpo-r   and  Zup&V  are  removed  as 
shown  in  Fig.  49,  Z^n  substitutes  for  Z-rz,  ZR^  substitutes 
for  Z-rp*  Ztt  substitutes  for  Z„,,  and  Z-ri,  substitutes  for 
Zrrp.   This  process  slurs  the  effect  that  would  be  present 
if  the  two  shunt  arms  were  not  removed  and  thus  improves 
the  complementary  filter  characteristics.   This  process 
will  only  be  applied  to  the  complementary  Zobel  filter  con- 
figuration. 


"LI 


ZL4  " 


'H4 


'L2 


;L3  J 


&H1 


'H2 


'H3 


Fig.  48.   Complementary  filter  configuration. 
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1XL 


Fig.  4-9.   Complementary  degenerate  filter  configuration 
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APPENDIX  B 


Aidentity  Algorithm.   Using  the  principle  of  King's  (2) 
approximate  identity  on  the  ADPI  of  a  filter  is  the  aidentity 
algorithm. 

Given  the  constant-1  filter  configuration  of  Pig.  50, 
the  aidentity  algorithm  is  applied  to  the  ADPI  of  equation 
(27).   This  process  is  to  eauate  as  many  successive  coeffi- 


ln 


2    5 

A  +A-,s+ApS  +A,s^ 

B  +B, s+B0s  +B2s^ 
old    p 


(27) 


where 


Ao  =  X 

Al  =  2 
A2  =  2  A 

A5  =  2  A  (2-  A) 


Bo  =  1 

Bl  =  2 


(28) 


B2  =  2(2-  X) 

B$  *  0 


cients  of  the  corresponding  powers  of  s  in  the  numerator 

and  denominator  as  is  possible.   This  has  A  =B  .  A,=B,  , 

o  o T  l  l  * 

and  Ap=Bp.   The  result  of  Ap  equated  to  Bp  is  Asl»   When 
this  value  of  A  is  substituted  in  A,=B^,  the  result  1^0 
implies  that  A,  and  B,  cannot  be  equated  if  Apand  Bp  are. 
Since  the  process  has  to  use  successive  coefficients,  \  =1 
will  be  used  and  the  resultant  ADPI  will  be  a  third  order 
ADPI.   However  this  ADPI  will  be  the  same  as  would  be  obtained 
if  the  constant-1  filter  of  Fig.  51  was  used.   Thus  this 
filter  aidentity  order  cannot  be  improved  by  the  aidentity 
algorithm. 


As 


zm  -t> 


(2-A)s 


-r   i/2s 


Fig.    50.      Low-pass   filter. 
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Pig.    51 •      Constant-1   filter. 
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APPENDIX  C 

Interactance.   Interactance  (1)  is  a  means  of  comparing 
the  power  available  at  the  input  terminals  of  a  filter, 
or  a  fan-out  filter  configuration  of  two  or  more  filters, 
with  that  which  would  be  available  if  all  of  the  filters 
were  purely  resistive.   Therefore  a  constant  interactance 
freouency  response  would  indicate  a  network  that  was  resistive 
in  nature.   The  ADPI  phase  shift  freouency  response  would 
give  the  same  information,  but  computational  wise  it  is  more 
difficult.   Interactance  is  more  realistic  than  character- 
istic impedance  in  fan-out  filter  configurations  because 
of  the  full  range  of  freouencies  that  are  being  encountered 
by  this  type  of  network.   Characteristic  impedance  only 
has  meaning  at  a  single  frequency. 

In  the  calculation  of  interactance, Vu  =Pf/Pr*  Pf  an<i 
P  are  calculated  for  e  being  a  peak  value  of  the  sinusoidal 
voltage  source  of  Pig.  52,  E  »«=|e|  Ay2,  and  the  source 
resistance,  R  ,  is  not  matched  to  the  filter  driving  point 
impedance.   Referring  to  Pig.  52,  P  is  the  power  dissipated 
in  the  source  resistance  when  a  pure  resistance  of  one  ohm 
replaces  each  filter  in  the  fan-out  network  and  a  voltage 
source,  -e/n,  is  placed  in  series  with  each  of  these  resist- 
ances.  Then  from  Fig.  52  using  Millman's  theorem,  R  =R1= 
.  .  .  =Rn=lA,  T0-l/R0,  VIg-  .  .  .  *Yn=lU< 
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e(l)  +  [(-e/n)(l)]n 


or 


1  +  (Dn 


(29) 


e   =  0 
or 


(30) 


e   =  e  -  e   =  e 
rr        or 


(3D 


Erreff  =  hrrl^" 


(32) 


Pr  =  Krl  /2Ro 


(33) 


.'.  Pr  =  h|2/2 


(3<0 


Referring  to  Fig.  53 »  ?*•  is  the  power  dissipated  in  the 

source  resistance  when  a  voltage  source,  -e/n,  is  placed 

in  series  with  each  fan-out  filter.   Then  from  Fig.  53  using 

Millman's  theorem,  R  =111 ,  Y  =1/R  . 

7   o       o     o 


of 


eYo  +  (-e/n)^  +  Y?  +  .  .  .  +  ?J 


Y^  +  Yn  +  Y~  +  .  .  .  +  Y^ 
o    1    c.  n 


(35) 


Tf 


e  -  eQf  =  e 


(n+l)(Y1+Y2+  .  .  .  +Yn) 


n(Yo+Y1+  .  .  .  +Yn)  _ 


(36) 


E 


Tf 


rfeff 


^ 


(57) 
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|2yo 


(38) 


.  •  -l  «  ■ 


(n+l)(y1+Y2+...+Yn) 


n(l+Y1+Y2+.,.+Yn) 


(39) 


The  interactance  is  now  the  ratio  of  P-  to  P  as  given  in 
eouation  (40). 


A- 


(n+l)(Y1+Y2+...+yn) 


n(l+Y1+Y2+ 


•*V 


(40) 


The  following  observations  should  be  made  concerning  the 
computational  technioues  just  given: 

1.  e,  e ,  e„.r>  *~.«i  and  ©~*  are  peak  values  of  sin- 

t   rr7   ri 7   or       of     F 

usoidel  voltages  and  used  for  ease  of  computation. 

2.  (-e/n)  voltages  were  added  for  ease  of  computation. 

3.  One  ohm  resistors  and  impedances  were  assumed  for 
ease  of  computation  since  filters  were  normalized 
to  impedance  level  of  one  ohm. 

4.  The  interactance  for  filters  at  different  impedance 
levels  can  be  determined,  but  this  just  confuses 
the  issue. 

As  an  example  to  show  the  interactance  calculations, 
the  constant-1  filter  is  used  as  shown  in  Fig.  51 •   The 
ADPI  is  given  in  equation  (41).   In  this  example  n=l  as 
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there  is  only  one  filter  end  Y,  is  the  inverse  of  the  ADPI 
as  shown  in  equation  (42).   Because  this  filter  is  normalized 
to  one  ohm,  Y  will  be  equal  to  one  mho.   Using  these  values 
in  eauation  (40)  results  in  equation  (43)  end   this  inter- 
actance  response  curve  is  given  in  Fig.  20. 


Z. 
in 


l+2s+2s2+2s5 


l+2s+2s' 


(<u) 


Yi  ■ 


l+2s+2s' 


l+2s+?s2+2s5 


(42) 


Since   Y  =1LT   and   n=l,    then 


A 


l+2s+2s' 


l+2s+2s2+s^ 


s=gw 


l+4w 


1+w' 


(43) 
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err>V1/flo 


or 


Fig.  52.   The  fan-out  resistive  network  used  for  the  deter- 
mination of  P  for  interactance. 


erf>V1/Ro 


Fig.  53.   The  fan-out  filter  network  used  for  the  determin- 
ation of  P«  for  interactance. 
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APPENDIX  D 

FORGO  Programs  for  Intcractencc  Calculations.   To  make 
possible  the  computation  of  the  interactance  versus  frequency 
graph,  the  FORGO  digital  computer  program  shown  on  pages 
65  through  69  was  written  and  used.  The  input  end  output 
data  for  each  of  the  previously  discussed  filter  configur- 
ations are  given  on  the  pages  following  the  computer  program. 


65 


READ  A(I) 
READ  B(I) 


I 


K=0 


I 


(1=10,120 


I 


H=0.0 
W=H/100, 
K=K+1 
W1(K)=V/ 


A 


=R(K) 


J=13,25; 


I 


P=J 

W=P/10. 
K=K+1 
V1(K)=W 


T 


A 


-R(K) 


4 


Q=M 

W=0/10 
K=K+1 
W1(K)=W 


Fig.  54.   Block  diagram  of  interactance  computer  program, 
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C   C   INTERACTANCE  VS.  FREQUENCY,  SINGLE  FILTERS 

DIMENSION  A(1G),B(10)»W1(140)»R(140) 
3  READ,AO, ( A(  I  )  ,  I  =  1  » 10  ) 

READ»Eo,(B(N) ,N=1,10 > 

K  =  0 

DC  2  1=10,120 

H=I 

W=H/100. 

K  =  K  +  1 

Wl  (K)=W 

D=(BO-B( 2)*W**2+B(4)*W**4-B(6)*W**6+B(8 )*W**8-B( 10 ) *W** 10 ) **2 

E=(B( 1 )*W-B( 3 )*W**3+B( 5 )*W**5-B(7) *W**7+B(9 )*W**9)**2 

Fl=A0+B0-( A( 2)+B( 2 ) )*W**2+( A (4 )+B( 4) )*W**4-( A( 6)+B(6) ) *W**6 

F2=( A(8 )+B( 8) )*W**8-( A( 10)+B( 10) )*W**10 

F= (F1+F2) **2 

Gl=(A(l)+B(l) )*W-( A( 3)+B(3 ) )*W**3+(A( 5 ) +B ( 5 ) ) *W**5 

G2=-(A(7)+B(7) )*W**7+(A(9)+B(9) )*W**9 

G= (G1+G2) **2 

R(K)=4.      *(D+E)/(F+G) 
2  CCNTINUE 

DC  5  J=13»25 

P=J 

W=P/10. 

K  =  K+1 

W1(K)=W 

D=(B0-B( 2)*W**2+B(4)*W**4-B(6 )*W**6+B(8 )*W**8-B( 10 ) *W**10 ) **2 

E=(B(1)*W-B(3 )*W**3+B( 5 )*W**5-B (7 ) *W**7+B ( 9 ) *W**9 ) **2 

Fl=A0+B0-( A( 2)+B(2 ) )*W**2+ ( A(4 )+B( 4) )*W**4-( A( 6 ) +B ( 6 ) )*W**6 

F2=( A(8 )+B(8) )*W**8-(A( 10)+B( 10) )*W**10 

F=(F1+F2)**2 

Gl= (A( 1 )+B( 1 ) )*W-( A( 3)+B(3 ) ',  *W**3+(  A( 5 )+B( 5) )*W**5 

G2=-(A(7)+B(7))*W**7+(A(9)+B(9) )*W**9 

G=(G1+G2)**2 

R(K)=4.      *(D+E)/(F+G) 
5  CCNTINUE 

DC  10  M=30,100,5 

Q  =  M 

W=Q/10. 

K  =  K  +  1 

W1(K)=W 

D=(B0-B(  2)*W**2+B(4)*W**4-B(6)*W**6+B(8  )*W**8-B(  10 ) *W**10 ) **2 

E=(B( 1 )*W-B(3)*W**3+B(5)*W**5-B(7)*W**7+B(9)*W**9)**2 

F1=A0+BG- ( A( 2)+B(2) )*W**2+ ( A(4)+B(4) )*W**4-( A( 6)+B( 6) ) *W**6 

F2=(A(8)+B(8) )*W**8~(A( 10)+B( 10) )*W**10 

F=(F1+F2)**2 

G1=(A(1)+B(1) )*W-(A(3)+B(3 ) ) *W**3+ ( A ( 5 ) +B ( 5) )*W**5 

G2=-(A(7)+B(7) )*W**7+( A(9)+B(9) )*W**9 

G=(G1+G2 ) **  2 

R(K)=4.      *(D+E)/(F+G) 
10  CCNTINUE 

PUNCH  4 
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DC 
PUN 


1 

CON 

DC 

PUN 

2 

CON 

GC 

4 

FOR 

7 

FOR 

8 

FOR 

11  K=l,45 

CH  7»Wl(K)tR(K)»Wl (K+47 )  ,R ( K+4  7 ) , Wl ( K  +  94 ) ,R ( K  +  94  ) 
TINUE 

12  K=46»47 

CH  8.W1 (K) ,R(K ) >W1 (K  +  47 ) ,R(K+47) 

TINUE 

TC  3 

MAT(/8X»1HW»9X,1HR,2( lOX » 1HW *9X ♦ 1HR ) / / ) 

MAT(3( 5X.F5.2»3X»F8.5)  ) 

MAT<2( 5X,F5.2»3X,F8.5)  ) 


END 


68 


C   C   INTERACTANCE  VS.  FREQUENCY,  COMPLEMENTARY  FILTERS 


10 


DIMENSION  A(10),B(10>,W1(140)»R(140) 

READ,AO,( A( I ) ,1=1,10) 

READ,BO,(B(N) ,N=1 ,10) 

K  =  0 

DO  2  1=10,120 

H=I 

W=H/100. 

K=K  +  1 

W1(K)=W 

D=(B0-B( 2 )*W**2+B(4)*W**4-B(6 )*W**6+B(8 )*W**8-B( 10 ) *W** 10 ) **2 

E=(B( 1 )*W-B( 3)*W**3+B(5)*W**5-B(7)*W**7+B(9)*W**9)**2 

Fl=AO+B0-( A( 2)+B(2 ) )*W**2+( A(4 )+B( 4) )*w**4-( A( 6)+3( 6) )*W**6 

F2=(A(8)+B(8) )*W**8-(A( 10)+B( 10) )*W**10 

F=(F1+F2)**2 

Gl=(A(l)+B(l) )*W-< A( 3)+B(3 ) )*W**3+(A( 5 )+B( 5) )*W**5 

G2  =  -(A(7)+B(7)  )*W**7+(  M9)+B(9  )  )*W**9 

G=(G1+G2)**2 

R  (K)=(9./4.  )*(D+F)  /(  F+G) 

CONTINUE 

DO  5  J=13,25 

P  =  J 

W=P/10. 

K  =  K+1 

Wl(K)=W 

D=(BO-B( 2)*W**2+B(4)*W**4-B(6 )*W**6+B( 8 )*W**8-B( 10)*W**10 )**2 

E=( B( 1 )*W-B( 3)*W**3+B( 5 )*W**5-B (7 )*W**7+B(9 )*W**9)**2 

Fl=A0+B0-( A( 2)+B(2) ) *W**2+ ( M4)+B(4) )*W**4-( A( 6)+B(6) )*W**6 

F2=(A(8)+B'8) )*W**8-(A( 10)+B( 10) )*W**10 

F=(F1+F2)*<>2 

Gl=(Ml)+B(l) )*W-( A(3)+R(3) )*W**3+(A( 5)+B( 5) )*W**5 

G2=-(A(7)+B(7) )*W**7+(A(9)+P(9) )»W**9 

G=(G1+G2)**2 

R(K)=(9./4.)*(D+E)/(F+G) 

CONTINUE 

DO  10  M=30, 100,5 

0  =  M 

W=0/10. 

K  =  K  +  1 

Wl (K)=W 

D=(BO-B( 2  5*W**2+B(4) *W**4-B ( 6 ) *W**6+B ( 8 )*W**8-B( 10 > *W**10 ) **2 

E=(B( 1 )*W-8( 3)*W**3+E ( 5 )*W**5~B (7) *W**7+B(9 )*W**9)**2 

Fl=A0+B0-( A( 2)+B(2 ) ) *W**2+ ( A ( 4 ) +B ( 4 ) ) *W**4- ( A ( 6 ) +B ( 6 ) ) *W**6 

F2=(A(8)+B(8) )*W**8-(A(10)+B( 10) )*W**10 

F=(F1+F2)**2 

G1=(A(  1  )+B(  1)  J*W-( A(  3)+B'(3)  )*W**3+(A(5)+B(5)  )*W**5 

G2  =  -(A(7)+B(7)) *W**7+( A(9) +R ( 9  )  )*W**9 

G=(G1+G2) **2 

R (K )=(9./4. )*<D+E )/( F+G) 

CONTINUE 

PUNCH  4 
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DC  11  K=1.45 

PUNCH  7,Wl(K)»R(K)*Wl(K+47)»R(K+4  7),Wl(K+94)  »R(K+94) 

11  CCNTINUE 

DC  12  K=46»47 

PUNCH  8,W1(K),R(K),W1 (K  +  47  )  ,R(K+4  7) 

12  CCNTINUE 
GC  TC  3 

4  FCRMAK /8X»1HW>9X»1HR,2(10X»1HW»9X*1HR)//) 

7  FCRMAT(3< 5X ♦ F5. 2 »3X , F8 . 5 ) ) 

8  FCRMAT(2(IX»F5.2»3X»F8.5) ) 
END 
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INTERACTANCE  RESPONSE  INPUT  DATA,  ADPI  CCEFF.,  CCNSTANT-1  FILTER 
1.   2.   2.   2.   .0   .0   .0   ,0   ,0   .0   .0 
la   2,   2.   .C   ,0   .0   .0   ,0   .0   .0   .0 


INTERACTANCE  RESPONSE  INPUT  DATA,  ADPI  CCEFF.,  ZOBEL  FILTER,  M=0.6 
1.   3.84   5.76   11.63   6.32   9.1   1.95   2.07   .0   .0   .0 
1.   3.2   8.06   6.64   10.81   3.07   3.26   .0   .0   .0   .0 


INTERACTANCE  RESPONSE  INPUT  DATA,  ADPI  CCEFF.,  ZOBEL  FILTER»  M=0.707 
la   3.414   6.828   9.246   8.992   6.682   2.914   1.457   .0   .0   .0 
la   3.414   6.828   9.242   8.992   5.828   2.914   .0   .0   .0   .0 


INTERACTANCE  RESPONSE  INPUT  DATA,  ADPI  CCEFF.,  TYPE  1  FILTER 
1.   3.   5.   6.5   5.   3.   1.   .0   .0   .0   .0 
1.   4.   8.   9.   8.   4.   1.   .0   .0   .0   .0 


INTERACTANCE  RESPONSE  INPUT  DATA,  ADPI  CCEFF.,  TYPE  2  FlLTEP 
3.24   11.7   21.55   28.79   21.55   11.7   3.24   .0   .0   .0   .0 
3.24   11.7   22.8   25.42   22.8   11.7   3.24   .0   .0   .0   .0 


INTERACTANCE  RESPONSE  INPUT  DATA,  ADPI  CCEFF.,  TYPE  3  FILTER 

3.73246   13.9295«   26.85862   36.3263   26.85862   13.92954   3.73246   .0 

.0   .0   .0 

3.73246   13.92954   26.8614   29.859   26.8614   13.92954   3.73246   .C   .0 

.0   .0 


INTERACTANCE  RESPONSE  INPUT  DATA,  ADPI  CCEFF.,  TYPE  4  FILTER 
.9975   3.99   7.98   9.95   7.98   3.99   .9975   .0   .0   .0   .0 
.9975   3.99   7.98   9.95   7.98   3.99   .9975   .0   .0   .0   .0 


INTERACTANCE  RE!>PCNSE  INPUT  DATA,  ADPI  CCEFF.,  TYPE  6  FILTER 
3.23607  15.70820  43.12461  *3. 33938  119.37381  135.99185  119.37381 
83.33938  43.12461  15.70820  3.23607         3.23607  15.70820  43.12461 
81.48529  118.37383  131.8459  118.37383  81.48529  43.12461  15.70820  3.23607 


INTERACTANCE  RESPCNSE  INPUT  DATA*  ADPI  CCEFF.,  TYPE  7  FILTER 

.68140           3.63470  10.42514  20.45234  29.90444 

33.93411  29.904444  20.45234  10.42514  3.6347 

.68140 

.68140           3.63470  10.42514  20.45234  29.93704 

33.77306  29.93704  20.45234  10.42514  3.6347 

.68140 
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INTERACTANCE  VS.  FREQUENCY,  OUTPUT  DATA,  CCNSTANT-1  FILTER,  FIG. 19 
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INTERACTANCE  VS.  FREQUENCY,  OUTPUT  DATA,  ZCBEL  FILTER,  FIG. 22 
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C   C   INTERACTANCE  VS.  FREQUENCY,  OUTPUT  DATA,  ZCBEL  FILTER,  FIG. 25 
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INTERACTANCE  VS.  FREQUENCY,  OUTPUT  DATA,  TYPE  1  FILTER,  FIG. 28 
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INTERACTANCE  VS.  FREQUENCY,  OUTPUT  DATA,  TYPE  2  FILTER,  FIG. 31 
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INTERACTANCE  VS.  FREQUENCY*  OUTPUT  DATA,  TYPE  3  FILTER,  FIG. 34 
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.72820 

1.05 

.11537 

.12 

.56323 

.59 

.73102 

1.06 

.12895 

.13 

.56350 

.60 

.73309 

1.07 

.14434 

.14 

.56384 

.6] 

.73436 

1  .08 

.16132 

.15 

.56426 

.62 

.73471 

1.09 

.17965 

.16 

.56477 

.63 

.73408 

1.10 

.19909 

.17 

.56537 

.64 

.73237 

1.11 

.21943 

.18 

•566C9 

.65 

.72948 

1.12 

.24044 

.19 

.56692 

.66 

.72535 

1.13 

.26194 

.20 

.56789 

.67 

.71988 

1.14 

.28373 

.21 

.56900 

.68 

.71299 

1.15 

.30563 

.22 

.57027 

.69 

.70462 

1.16 

.32751 

.23 

.57170 

.70 

.69470 

1.17 

.34922 

.24 

.57331 

.71 

.68318 

1.18 

.37064 

.25 

.57511 

.72 

.67000 

1.19 

.39167 

.26 

.57710 

.73 

.65515 

1.20 

.41222 

.27 

.57930 

.74 

.63861 

1.30 

.58068 

.28 

.58171 

.75 

.62038 

1.40 

.67774 

.29 

.58435 

.76 

.60048 

1.50 

.72186 

.30 

.58721 

.77 

.57897 

1.60 

.73453 

.31 

.59031 

.78 

.55592 

1.70 

.73057 

.32 

.59365 

.79 

.53143 

1.80 

.71875 

.33 

.59724 

.80 

.50561 

1.90 

.70390 

.34 

.60107 

.81 

.47863 

2.00 

.68855 

.35 

.60514 

.82 

.45067 

2.10 

.67392 

.36 

.60946 

.83 

.42193 

2.20 

.66057 

.37 

.61401 

.84 

.39266 

2.30 

.64866 

.38 

.61880 

.85 

.36312 

2.40 

.63818 

.39 

.62381 

.86 

.33359 

2.50 

.62904 

.40 

.62904 

.87 

.30438 

3.00 

.59849 

.41 

.63446 

.88 

.27578 

3.50 

.58319 

.42 

.64007 

.89 

.24813 

4.00 

.57511 

.43 

.64585 

.90 

.22173 

4.50 

.57057 

.44 

.65176 

.91 

.19691 

5.00 

.56789 

.45 

.65780 

.92 

.17394 

5.50 

.56623 

.46 

.66392 

.93 

.15310 

6.00 

.56516 

.47 

.67009 

.94 

.13464 

6.50 

.56444 

.48 

.67628 

.95 

.11876 

7.00 

.56395 

.49 

.68245 

.96 

.10561 

7.50 

.56361 

.50 

.68855 

.97 

.09533 

8.00 

.56336 

.51 

.69453 

.98 

.08798 

8.50 

.56317 

.52 

.70034 

.99 

.08359 

9.00 

.56304 

.53 

.70593 

1.00 

.08214 

9.50 

.56293 

.54 

.71123 

1.01 

.08356 

10.00 

.56285 

.55 

.71618 

1.02 

.08775 

.56 

.72071 

1.03 

.09457 
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C   C   INTERACTANCE  VS.  FREQUENCY*  OUTPUT  DATA,  TYPE  6  FILTER,  FIG. 42 
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.58870 
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.40035 
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.58784 
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.43968 
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.56300 

.61 

.58684 
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.47297 
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.56315 
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.58572 

1.09 

.50020 
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.56333 

.63 

.58448 

1.10 

.52186 

.17 

.56355 
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.58315 

1.11 

.53865 

.18 

.56380 

.65 

.58174 

1.12 

.55134 

.19 

.56409 

.66 

.58028 

1.13 

.56068 

.20 

.56442 

.67 

.57881 

1.14 

.56734 

.21 

.56480 

.68 

.57734 

1.15 

.57190 

.22 

.56523 

.69 

.57592 

1.16 

.57483 

.23 

.56570 

.70 

.57459 

1.17 

.57655 
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.56623 

.71 

.57338 
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.57736 
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.56680 
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.57233 

1.19 

.57752 

.26 

.56743 

.73 

.57149 

1.20 

.57724 

.27 

.56812 

.74 

.57089 

1.30 

.57080 

.28 

.56885 

.75 

.57057 

1.40 

.57291 

.29 

.56964 

.76 

.57055 

1.50 

.57930 

.30 

.57048 

.77 

.57084 

1  .60 

.58511 

.31 

.57137 

.78 

.57144 

1.70 

•588P4 

.32 

.57230 

.79 

.57233 

1  .80 

.59052 

.33 

.57327 

.80 

.57345 

1.90 

.59067 

.34 

.57429 

.81 

.57471 

2.00 

.58979 

.35 

.57534 

.82 

.57596 

2.10 

.58831 

.36 

.57641 

.83 

.57699 

2.20 

.58652 

.37 

.57751 

.84 

.57752 

2.30 

.58460 

.38 

.57862 

.85 

.57716 

2.40 

.58269 

.39 

.57974 

.86 

.57542 

2.50 

.58086 

.40 

.58086 

.87 

.57168 

3.00 

.57361 

.41 

.58196 

.88 

.56519 

3.50 

.56930 

.42 

.58305 

.89 

.55504 

4.00 

.56680 

.43 

.58411 

.90 

.54024 

4.50 

.56533 

.44 

.58513 

.93 

.51973 

5.00 

.56442 

.45 

.58610 

.92 

.49253 

5.50 

.56385 

.46 

.58700 

.93 

.45799 

6.00 

.56347 

.47 

.58784 

.94 

.41609 

6.50 

.56322 

.48 

.58859 

.95 

.36797 

7.00 

.56304 

.49 

.58924 

.96 

.31638 

7.50 

.56291 

.50 

.58979 

.97 

.26597 

8.00 

.56282 

.51 

.59023 

.98 

.22296 

8.50 

•56275 

.52 

.59054 

.99 

.19391 

9.00 

.56270 

.53 

.59072 

1.00 

.18366 

9.50 

.56266 

.54 

.59075 

1.01 

.19369 

10.00 

.56263 

.55 

.59065 

1.02 

.22152 

.56 

.59039 

1.03 

.26180 
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INTERACTANCE  VS.  FREQUENCY,  OUTPUT  DATA,  TYPE  7  FILTER,  FIG. 45 
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.56226 

.75 

.62880 

1  .40 

.60448 

.29 

.56220 

.76 

.63641 

1.50 

.58114 

.30 

.56213 

.77 

.64422 

1  .60 
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.56131 

.84 
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.56007 

.38 

.56115 

.85 
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2.40 

.56045 

.39 

.56097 

.86 

.68351 
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.56078 

.40 

.56078 

.87 
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.56182 

.41 

.56058 
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.42 

.56038 
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.56017 
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4.50 

.56245 
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.55997 
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.46262 
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.56250 
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.55917 

.96 

.35318 
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.56250 

.50 
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.97 

.30224 
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.56250 

.51 

.55911 

.96 

.26067 

8-50 

.56250 
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.55919 

.99 

.23348 
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.56250 
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.55934 

1.00 
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9.50 

.56250 
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.55960 

1.01 
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10.00 

.56250 

.55 

.55998 

1.02 

.25931 

.56 

.56050 

1.03 

.29815 
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APPENDIX  E 

FORGO  Programs  for  Aidentity  Driving  Point  Impedance 
(ADPI)  Calculations.   To  make  possible  the  computation  of  the 
ADPI  versus  freauency  graph,  the  FORGO  digital  computer 
program  shown  on  pages  80  and  81  was  written  and  used. 
The  input  end  output  data  for  each  of  the  previously  dis- 
cussed filter  configurations  are  given  on  the  pages  following 
the  computer  program.   The  program  would  not  handle  Type  5* 
a  fourteenth  order  polynomial,  since  it  had  a  maximum  cap- 
acity for  a  tenth  order  polynomial. 


80 


r 


READ  A(I) 
READ   B(I) 


K=0 


I 


1=10,12) 


i 


K=K+1 
H=I 

W=H/100. 
tfl(K)-W 


I 


ZOO- 


Z(s) 


s=dwl(K) 


J. 


THE(K)=ARGZ(K) 


PUNCH 

COLUMN 

HEADINGS 


l 


K-1,11 


I 


PUNCH 

ZOO 

W1(K) 
THS(K) 


Fig.  55.   Block  diagram  of  ADPI  computer  program. 


81 


ADPI(  Z(S)  -  PHASE  SHIFT  )  VS.  FREQUENCY,  ALL  FILTERS 


4 

5 

30 

31 

45 


DIM 
FOR 
FOR 
FOR 
FOR 
REA 
K  =  0 
DC 
K  =  K 
H=I 
W  =  H 
Wl  ( 
C  =  A 
D  =  A 
E  =  B 
F  =  B 
Z(K 
IF( 
3  PUN 
THE 
GO 
8  I  F( 
11  IF( 
13  THE 


10 
16 

15 

14 

9 


22 

19 
25 

24 

20 

2 


ENSIGN  A( 
MAT( /7X,2 
MAT(  2(5X 
MAT( 7HW(K 
MAT( 7HW(K 
D,AO, ( A(  I 

2  1=10,12 
+  1 

/10-. 
K)=W 

0-A(2 ) *W* 
( 1 )*W-A( 3 
0-B(2 ) *W* 
( 1 )*to-B(3 
)=SQRT( (C 
ABS(C)-.l 
CH  30,W1( 
1  =  0. 
TC  14 
C  )  1 0  , 1  0  , 1 
D)12,12,l 
1=ATAN(D/ 


10),B(10),W1(112),Z(112),THE(112) 
!HW1>5X,5HZ(JW),6X»2HC1»9X,2HW1,5X,5HZ(JW),6X,2HC1// 

,F5.2,3X,F7.5,2X,F7.2)  ) 

.)     =    F5.2,5X,4HC    =    E15.8,5x,4HD    =    E15.8) 
.)     =    F5.2,5X,4HE    =    El5.8,5x,4HF    =    E15.8) 

) , 1=1,10) ,B0, ( B( I ) , 1=1,1^) 


)     =    F5.2,5X,4HE    =    El5.o,5X 
), 1=1,10) , BO, ( B( I ), 1=1,10) 

0 


*2+A(4)*W**4-A(6)*W**6+A(8 )*W**8-A< 10 )*W**10 
)*W**3+A ( 5)*W**5-A(7)*W**7+A( 9)*w**9 
*2+B(4)*W**4-B ( 6 )*W**6+B( 8 )*W**8-B( 10 )*W**10 
)*W**3+B(  5)*W**5-B(7)*W**7+B(  9)*W**9 
**2+D**2 )/(E**2+F**2) ) 
E-40)3,3,8 
K)  ,C,D 


0*57.245779 
0/0*57.24779 


GO    TO    14 


GO    TO    14 


IF( 
THE 
GO 
THE 
IF( 
PUN 
THE 
GO 
28     IF( 
21     IF( 
23    THE 
GO 


6 


THE 

GO 

IF( 

THE 

GO 

THE 

THE 

CON 

PUN 

DO 

PUN 

CON 


D)15,15,li 

1=180.-ATAN(D/ABS(C))*5  7.245779 
TO    14 

1  =  180.  +AT- 
ABS(E)-.l 
CH    31,W1( 

2  =  0. 
TO    20 
E)19,19,2 
F)22,22,2: 
2=ATAN(F/I 
TO    20 
2=ATAN (F/l 
TO    2  0 

F)24,24,25 

2=180.-ATAN(F/ABS(E) )*5  7.2  45  779 
TO  20 
2=180. +AT, 
<K)=THE1-- 
TINUE 

CH  4 

7  J=l,55 

CH  5,Wl{J),Z(J),THE(J),Wl(J+55)  ,Z(J+55),THE(J+55) 

TINUE 


AN(D/C)*57. 245779 

E-40)9,9,28 

K) ,E,F 


1 
3 
E)*57. 245779 

E)*57. 245779 


AN(F/E)*57. 245779 
THE2 


GO  T: 
END 
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ADPI  RESPONSE,  ADPI  COEFFICIENTS  -  INPUT  DATA,  CONSTANT-1  FILTER 
2.   2.   2.   .0   .0   .0   .0   .0   .0   .0 

2.  2.   .0   .0   .0   .0   .0   .0   .0   .0 

ADPI  RESPONSE,  ADPI  COEFFICIENTS  -  INPUT  DATA,  ZOBEL  FILTER,  M=0.6 
3.84   5.76   11.63   6.32   9.1   1.95   2.07   .0   .0   .0 
3.2   8.06   6.64   10.81   3.07   3.26   .0   .0   .0   .0 

ADPI  RESPONSE,  ADPI  COEFFICIENTS  -  INPUT  DATA,  ZOBEL  FILTER,  M=0.707 
3.414   6.828   9.246   8.992   6.682   2.914   1.457   .0   .0   .0 
3.414   6.828   9.242   8.992   5.828   2.914   .0   .0   .0   .0 

ADPI  RESPONSE,  ADpI  COEFFICIENTS  -  INPUT  DATA,  TYPE  1  FILTER 

3.  5.   6.5   5.   3.   1.   .0   .0   .0   .0 

4.  8.   9.   8.   4.   1.   .0   .0   .0   .0 

ADPI  RESPONSE,  ADPI  COEFFICIENTS  -  INPUT  DATA,  TYPE  2  FILTER 
3.24   11.7   21.55   28.79   21.55   11.7   3.24   .0   .0   .0   .0 
3.24   11.7   22.8   25.42   22.8   11.7   3.24   .0   .0   .0   .0 

ADPI  RESPONSE,  ADPI  COEFFICIENTS  -  INPUT  DATA,  TYPE  3  FILTER 

3.73246   13.92954   26.85862   36.3263   26.85862   13.92954   3.73246   .0 

•  0   .0   •  C 

3.73246   13.92954   26.8614   29.859   26.8614   13.92954   3.73246   .0   .0 

•  0   .0 

ADPI  RESPONSE,  ADPI  COEFFICIENTS  -  INPUT  DATA,  TYPE  4  FILTER 
.9975   3.99   7.98   9.95   7.98   3.99   .9975   .0   .0   .0   .0 
.9975   3.99   7.98   9.95   7.98   3.99   .9975   .0   .0   .0   .0 

ADPI  RESPONSE,  ADPI  COEFFICIENTS  -  INPUT  DATA,  TYPE  6  FILTER 
3.23607  15.70820  43.12461  83.33938  119.37381  135.99185  119.37381 
83.33938  43.12461  15.70820  3.23607         3.23607  15.70820  43.12461 
81.48529  118.37383  131.8459  118.37383  81.48529  43.12461  15.70820  3.23607 

ADPI  RESPONSE,  ADPI  COEFFICIENTS  -  INPUT  DATA,  TYPE  7  FILTER 

.68140  3.63470  10.42514  20.45234  29.90444 

33.93411  29.904444  20.45234  10.42514  3.6347 
.68140 

.68140  3.63470  10.42514  20.45234  29.93704 

33.77306  29.93704  20.45234  10.42514  3.6347 
•6814C 


83 
C   ADPI (Z(S)-PHASE)  VS.  FREQUENCY,  OUTPUT  DATA,  CCNSTANT-1  FILTER,  FIG. 19 
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Z(  JW) 


CI 


Wl 


Z  (  JW  ) 


CI 


.10 

.99960 

-.11 

.65 

.58553 

-4.86 

.11 

.99942 

-.15 

.66 

.57007 

-4.23 

.13 

.99887 

-.24 

.68 

.53961 

-2.7] 

.14 

.99848 

-.30 

.69 

.52473 

-1.80 

.15 

.99800 

-.37 

.70 

.51015 

-.79 

.16 

.99742 

-.44 

.71 

.49593 

.34 

.17 

.99671 

-.53 

.72 

.48212 

1.57 

.18 

.99588 

-.62 

.73 

.46879 

2.93 

.19 

.99489 

-.73 

.74 

.45601 

4.41 

.21 

.99242 

-.97 

.76 

.43232 

7.75 

.22 

.99090 

-1.10 

.77 

.42154 

9.62 

.23 

.98916 

-1.25 

.78 

.41157 

11.62 

.24 

.98720 

-1.40 

.79 

.40247 

13.74 

.25 

.98498 

-1.56 

.80 

.39429 

15.98 

.26 

.98250 

-1.74 

.81 

.38709 

18.34 

.27 

.97974 

-1.92 

.82 

.38093 

20.80 

.28 

.97667 

-2.12 

.83 

.37585 

23.35 

.29 

.97329 

-2.32 

.84 

.37188 

25.96 

.30 

.96957 

-2.53 

.85 

.36905 

28.63 

.31 

.96549 

-2.75 

.86 

.36736 

31.34 

.32 

.96105 

-2.98 

.87 

.36682 

34.06 

.33 

.95622 

-3.21 

.88 

.36742 

36.76 

.35 

.94535 

-3.70 

.90 

.37193 

42.08 

.36 

.93928 

-3.95 

.91 

.37575 

44.64 

.37 

.93278 

-4.20 

.92 

.38056 

47.13 

.38 

.92582 

-4.46 

.93 

.38630 

49.53 

.39 

.91841 

-4.71 

.94 

.39290 

51.84 

.40 

.91054 

-4.97 

.95 

.40031 

54.03 

.41 

.90220 

-5.22 

.96 

.40846 

56.12 

.42 

.89339 

-5.47 

.97 

.41728 

58.10 

.44 

.87435 

-5.95 

.99 

.43671 

61.73 

.45 

.86412 

-6.18 

1.00 

.44721 

63.38 

.46 

.85344 

-6.40 

1.01 

.45816 

64.93 

.47 

.84229 

-6.61 

1.02 

.46952 

66.39 

.48 

.83070 

-6.80 

1.03 

.48123 

67.75 

.49 

.81867 

-6.97 

1.04 

.49326 

69.03 

.50 

.80623 

-7.12 

1.05 

.50557 

70.22 

.51 

.79338 

-7.25 

1.06 

.51813 

71.33 

.52 

.78014 

-7.35 

1.07 

.53091 

72.38 

.53 

.76654 

-7.43 

1.08 

.54388 

73.35 

.54 

.75260 

-7.47 

1.09 

.55701 

74.26 

.55 

.73835 

-7.48 

1.10 

.57028 

75.12 

.56 

.72381 

-7.45 

1.11 

.58368 

75.91 

.58 

.69398 

-7.27 

1.13 

.61076 

77.36 

.59 

.67876 

-7.11 

1.14 

.62442 

78.02 

.60 

.66339 

-6.89 

1.15 

.63815 

78.63 

.61 

.64789 

-6.62 

1.16 

.65191 

79.21 

.62 

.63232 

-6.28 

1.17 

.66572 

79.75 

.63 

.61671 

-5.88 

1.18 

.67956 

80.25 

.64 

.60109 

-5.41 

1.19 

.69341 

80.73 
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C   C   ADPI (Z(S)-PHASE)  VS.  FREQUENCY,  OUTPUT  DATA,  ZCBEL  EILTER,  FIG. 22 


Wl 


Z(  JW) 


Wl 


Z(  JW) 


01 


.10    ] 

.04274 

2.75 

.65 

.54949 

10.78 

.12    ] 

..05972 

2.84 

.67 

.56233 

17.34 

.13    ] 

.06887 

2.80 

.68 

.57453 

20.59 

.15    ] 

[  .08798 

2.54 

.70 

.61145 

26.77 

.16    ] 

..09774 

2.32 

.71 

.63658 

29.63 

.18    ] 

[.11705 

1.68 

.73 

.70144 

34.75 

.19    ] 

L. 12634 

1.25 

.74 

.74188 

36.97 

.20    ] 

..13519 

.77 

.75 

.78831 

38.93 

.21    1 

..14346 

.21 

.76 

.84138 

40.64 

.22    ] 

..15098 

-.40 

.77 

.90198 

42.07 

.23    ] 

..15761 

-1.08 

.78 

.97125 

43.23 

.25    ] 

L. 16755 

-2.60 

.80 

1.14217 

44.63 

.26    ] 

..17055 

-3.44 

.81 

1.24835 

44.82 

.27    1 

L.17206 

-4.32 

.82 

1.37262 

44.61 

.29    ] 

L. 17011 

-6.20 

.84 

1.69520 

42.68 

.30    ] 

..16644 

-7.18 

.85 

1.90766 

40.71 

.31    1 

[.16089 

-8.18 

.86 

2.16725 

37.80 

.32    ] 

L. 15341 

-9.18 

.87 

2.48585 

33.60 

.33    ] 

[.14398 

-10.18 

.88 

2.87282 

27.63 

.34    ] 

L. 13262 

-1  1.18 

.89 

3.32127 

1°.23 

.35    ] 

L. 11938 

-12.16 

.90 

3.77532 

7.77 

.36    : 

[.10430 

-13.11 

.91 

4.09281 

-6.76 

.37    ] 

[.08749 

-14.02 

.92 

4.10160 

-22.88 

.38 

L. 06905 

-14.88 

.93 

3.78655 

321.88 

.39 

L. 04911 

-15.69 

.94 

3.31218 

309.64 

.40    ] 

L. 02783 

-16.59 

.95 

2.83484 

300.43 

.41 

L. 00535 

-17.27 

.96 

2.42098 

293.66 

.42 

.98185 

-17.85 

.97 

2.08054 

288.67 

.43 

.95748 

-18.35 

.98 

1.80385 

284.92 

.44 

.93244 

-18.75 

.99 

1.57788 

282.05 

.45 

.90688 

-19.04 

1.00 

1.39124 

279.80 

.46 

.88098 

-19.21 

1.01 

1.23501 

278.02 

.47 

.85491 

-19.25 

1.02 

1.10248 

276.58 

.48 

.82884 

-19.00 

1.03 

.98861 

275.40 

.49 

.80292 

-18.77 

1.04 

.88962 

274.42 

.50 

.77731 

-18.38 

1.05 

.80264 

273.60 

.51 

.75218 

-17.83 

1.06 

.72546 

272.89 

.52 

.72769 

-17.10 

1.07 

.65638 

272.28 

.53 

.70398 

-16.19 

1.08 

.59406 

271.75 

.54 

.68123 

-15.09 

1.09 

.53743 

271.28 

.55 

.65961 

-13.78 

1.10 

.48562 

270.86 

.56 

.63929 

-12.26 

1  .11 

.43796 

270.48 

.57 

.62046 

-10.53 

1.]? 

.39386 

270.12 

.58 

.60333 

-8.56 

1.13 

.35286 

26°. 7Q 

.59 

.58809 

-6.38 

1.14 

.31457 

269.47 

.60 

.57498 

-3.97 

1.15 

.27866 

269.15 

.61 

.56422 

-1.36 

1.16 

.24486 

268.82 

.62 

.55605 

1.46 

1.17 

.21292 

268.49 

.63 

.55072 

4.44 

1.18 

.18265 

268.12 

.64 

.54845 

7.56 

1.19 

.15388 

267.69 

85 
C   C   ADPI (Z(S)-PHASE)  VS.  FREQUENCY,  OUTPUT  DATA,  ZCBEL  FILTER,  FIG. 25 


Wl 


Z<  JW) 


CI 


Wl 


Z  (  JW  ) 


.10 

1.00000 

.00 

.65 

1.02577 

-1.96 

.12 

1.00001 

.00 

.67 

1.02217 

-2.11 

.13 

1.C0001 

.00 

.68 

1.01987 

-2.14 

.15 

1.00002 

.00 

.70 

1.01446 

-2.07 

.17 

1.00005 

.00 

.72 

1.00851 

-1.77 

.19 

1.00011 

.00 

.74 

1.00312 

-1.16 

.20 

1.00015 

.01 

.75 

1.00116 

-.72 

.22 

1.00026 

.02 

.77 

1.00028 

.52 

.23 

1.00034 

.02 

.78 

1.00221 

1.33 

.24 

1.00044 

.02 

.79 

1.00644 

2.27 

.25 

1.000  5  5 

.03 

.80 

1.01368 

3.34 

.26 

1.00070 

.03 

.81 

1.02477 

4.56 

.27 

1.00086 

.04 

.82 

1.04076 

5.89 

.28 

1.00106 

.04 

.83 

1.06290 

7.34 

.29 

1.00130 

.05 

.84 

1.09277 

8.86 

.30 

]  .00157 

.05 

.85 

1.13232 

10.40 

.31 

1  .00188 

.06 

.86 

1.18399 

11.92 

.32 

1.00224 

.06 

.87 

1.25089 

13.32 

.33 

1.00264 

.07 

.88 

1.33705 

14.50 

.34 

1.00  3  10 

.07 

.89 

1.44774 

15.32 

.35 

1.00362 

.07 

.90 

1.58996 

15.58 

.36 

1.00420 

.07 

.91 

1.77286 

15.03 

.37 

1.00484 

.07 

.92 

2.00788 

13.34 

.38 

1.00554 

.07 

.93 

2.30713 

9.99 

.39 

1.00631 

.07 

.94 

2.67559 

4.33 

.40 

1.00716 

.06 

.95 

3.08795 

-4.45 

.41 

1.00807 

.05 

.96 

3.44631 

-16.76 

.42 

1.00905 

.04 

.97 

3.58132 

-31.62 

.43 

1.01010 

.02 

.98 

3.40664 

-46.37 

.44 

1.01121 

-.00 

.99 

3.03443 

-58.57 

.45 

1.01238 

-.03 

1.00 

2.62587 

-67.56 

.46 

1.01361 

-.07 

1.01 

2.26314 

-73.89 

.47 

1.01489 

-.11 

1.02 

1.96429 

-78.30 

.48 

1.01621 

-.16 

1.03 

1.72282 

-81.42 

.49 

1.01755 

-.21 

1.04 

1.52714 

-83.64 

.50 

1.0189] 

-.28 

1.05 

1.36673 

-85.24 

.51 

1.02027 

-.35 

1.06 

1.23332 

-86.43 

.52 

1.02161 

-.43 

1.07 

1.12073 

-87.30 

.53 

1.02291 

-.52 

1.08 

1.02435 

-87.96 

.54 

1.02415 

-.61 

1.09 

.94078 

-88.46 

.55 

1.02530 

-.72 

1.10 

.86745 

-88.84 

.56 

1.02634 

-.83 

1.11 

.80243 

-89.13 

.57 

1.02723 

-.95 

1.12 

.74422 

-89.36 

.58 

1.02796 

-1.08 

1.13 

.69165 

-89.54 

.59 

1.02848 

-1.21 

1.14 

.64381 

-89.67 

.60 

1.02877 

-1.34 

1.15 

.59996 

-89.78 

.61 

1.02880 

-1.48 

1.16 

.55954 

-89.87 

.62 

1.02854 

-1.61 

1.17 

.52205 

-89.93 

.63 

1.02796 

-1  .74 

1.18 

.48711 

-89.99 

.64 

1.02704 

-1.86 

1.19 

.45439 

-90.03 

86 
C   C   ADPI (Z(S)-PHASE)  VS.  FRFQUFNCY,  OUTPUT  DATA,  TYPF  1  FILTER,  FIG. 28 


Wl 


Z(  JW) 


Wl 


2(  JW) 


CI 


.10 

.99440 

-5.84 

.65 

.49916 

-30.07 

.12 

.99157 

-7.07 

.67 

.47600 

-28.91 

.15 

.98581 

-8.96 

.70 

.44404 

-26.53 

.18 

.97785 

-10.92 

.73 

.41636 

-23.32 

.19 

.97463 

-11.59 

.74 

.40831 

-22.05 

.20 

.97111 

-12.27 

.75 

.40093 

-20.69 

.21 

.96726 

-12.95 

.76 

.39430 

-19.24 

.22 

.96307 

-13.64 

.77 

.38847 

-17.69 

.23 

.95853 

-14.33 

.78 

.38351 

-16.06 

.24 

.95362 

-15.03 

.79 

.37950 

-14.36 

.25 

.94834 

-15.73 

.80 

.37651 

-12.59 

.26 

.94266 

-16.43 

.81 

.37462 

-10.77 

.27 

.93659 

-17.14 

.82 

.37389 

-8.93 

.28 

.93012 

-17.85 

.83 

.37440 

-7.09 

.29 

.92323 

-18.55 

.84 

.37620 

-5.26 

.30 

.91594 

-19.26 

.85 

.37933 

-3.49 

.31 

.90822 

-19.96 

.86 

.38383 

-1.81 

.32 

.90010 

-20.66 

.87 

.38968 

-.24 

.33 

.89156 

-21.35 

.88 

.39687 

1.18 

.34 

.88262 

-22.04 

.89 

.40530 

2.41 

.35 

.87328 

-22.72 

.90 

.41487 

3.43 

.36 

.86355 

-23.38 

.91 

.42537 

4.21 

.37 

.85344 

-24.04 

.92 

.43655 

4.71 

.38 

.84297 

-24.68 

.93 

.44808 

4.95 

.39 

.83215 

-25.46 

.94 

.45953 

4.90 

.40 

.82099 

-26.07 

.95 

.47044 

4.57 

.41 

.80952 

-26.66 

.96 

.48029 

3.99 

.42 

.79776 

-27.23 

.97 

.48857 

3.20 

.43 

.78571 

-27.78 

.98 

.49482 

2.23 

.44 

.77342 

-28.31 

.99 

.49869 

1.14 

.45 

.76089 

-28.80 

1.00 

.50000 

0.00 

.46 

.74816 

-29.27 

1.01 

.49872 

-1.13 

.47 

.73524 

-29.72 

1.02 

.49502 

-2.19 

.48 

.72216 

-30.12 

1.03 

.48920 

-3.12 

.49 

.70895 

-30.50 

1.04 

.48168 

-3.89 

.50 

.69563 

-30.84 

1.05 

.47290 

-4.46 

.51 

.68222 

-31.14 

1.06 

.46332 

-4.82 

.52 

.66875 

-31.25 

1.07 

.45336 

-4.96 

.53 

.65524 

-31.47 

1.08 

.44337 

-4.89 

.54 

.64172 

-31.65 

1.09 

.43363 

-4.61 

.55 

.62821 

-31.79 

1.10 

.42439 

-4.15 

.56 

.61474 

-31.87 

1.11 

.41578 

-3.51 

.57 

.60133 

-31 .91 

1.12 

.40793 

-2.73 

.58 

.58800 

-31.89 

1.13 

.40090 

-1.82 

.59 

.57478 

-31.82 

1.14 

.39472 

-.80 

.60 

.56170 

-31.69 

1.15 

.38940 

.30 

.61 

.54877 

-31.50 

1.16 

.38493 

1.47 

.62 

.53603 

-31.25 

1.17 

.38127 

2.69 

.63 

.52349 

-30.92 

1.18 

.37841 

3.94 

.64 

.51119 

-30.53 

1.19 

.37628 

5.20 
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C   C   ADPI (Z(S)-PHASE)  VS.  FREQUENCY,  OUTPUT  DATA,  TYPE  2  FILTER,  FIG. 31 
Wl      Z(JW)       01  Wl      Z(JW)       01 


.10 

1.00322 

-.13 

.65 

.76985 

12.63 

.13 

1.00471 

-.28 

.68 

.78165 

17.25 

.15 

1.00550 

-.43 

.70 

.79891 

20.48 

.17 

1.00595 

-.60 

.72 

.82508 

23.71 

.19 

1.00589 

-.81 

.74 

.86144 

26.85 

.20 

1.00562 

-.93 

.75 

.88389 

28.34 

.21 

1.00516 

-1.06 

.76 

.90944 

29.77 

.22 

1.00448 

-1.19 

.77 

.93830 

31.12 

.23 

1.00357 

-1.33 

.78 

.97070 

32.37 

.24 

1.00241 

-1.48 

.79 

1.00690 

33.52 

.25 

1.00097 

-1.63 

.80 

1.04717 

34.54 

.26 

.99923 

-1.78 

.81 

1.09181 

35.43 

.27 

.99718 

-1.94 

.82 

1.14118 

36.16 

.28 

.99480 

-2.^9 

.83 

1.19562 

36.73 

.29 

.99206 

-2.25 

.84 

1.25554 

37.1] 

.30 

.98897 

-2.41 

.85 

1.32136 

37.29 

.31 

.98550 

-2.56 

.86 

1.39350 

37.25 

.32 

.98164 

-2.71 

.87 

1.47237 

36.96 

.33 

.97739 

-2.85 

.88 

1.55832 

36.40 

.34 

.97274 

-2.98 

.89 

1.65157 

35.54 

.35 

.96769 

-3.10 

.90 

1.75214 

34.34 

.36 

.96225 

-3.21 

.91 

1.85966 

32.78 

.37 

.95640 

-3.30 

.92 

1.97323 

30.82 

.38 

.95017 

-3.37 

.93 

2.09113 

28.43 

.39 

.94357 

-3.42 

.94 

2.21061 

25.58 

.40 

.93660 

-3.44 

.95 

2.32763 

22.25 

.41 

.92929 

-3.44 

.96 

2.43681 

18.47 

.42 

.92165 

-3.56 

.97 

2.53171 

14.26 

.43 

.91372 

-3.34 

.98 

2.60555 

9.70 

.44 

.90552 

-3.23 

.99 

2.65236 

4.90 

.45 

.89710 

-3.09 

1.00 

2.66832 

0.00 

.46 

.88847 

-2.90 

1.01 

2.65267 

-4.35 

.47 

.87970 

-2.66 

1.02 

2.60793 

-9.52 

.48 

.87082 

-2.38 

1.03 

2.53910 

-13.88 

.49 

.86189 

-2.04 

1.04 

2.45251 

-17.85 

.50 

.85295 

-1.64 

1.05 

2.35456 

-21 .40 

.51 

.84407 

-1.19 

1  .06 

2.25087 

-24.50 

.52 

.83532 

-.67 

1.07 

2.14586 

-27.18 

.53 

.82675 

-.08 

1.08 

2.04272 

-29.46 

.54 

.81844 

.57 

1.09 

1.94355 

-31.37 

.55 

.81048 

1.30 

1.10 

1.84962 

-32.94 

.56 

.80293 

2.09 

1.11 

1.76155 

-34.22 

.57 

.79589 

2.97 

1.12 

1.67957 

-35.23 

.58 

.78945 

3.92 

1.13 

1.60361 

-36.01 

.59 

.78370 

4.94 

1.14 

1.53346 

-36.59 

.60 

.77875 

6.05 

1.15 

1.46879 

-36.98 

.61 

.77470 

7.22 

1.16 

1.40925 

-37.21 

.62 

.77167 

8.47 

1.17 

1.35446 

-37.30 

.63 

.76977 

9.80 

1.18 

1.30404 

-37.27 

.64 

.76912 

11.18 

1.19 

1.25766 

-37.12 

88 
C   ADPI (Z(S)-PHASE)  VS.  FREQUENCY,  OUTPUT  DATA,  TYPE  3  FILTER,  FIG. 34 


Wl 


Z(  JW) 


01 


Wl 


Z(  JW) 


01 


.10 

.99937 

-.09 

.65 

.80153 

25.89 

.13 

.99823 

-.19 

.68 

.84372 

31.14 

.15 

.99689 

-.29 

.70 

.88441 

34.53 

.17 

.99493 

-.39 

.72 

.93662 

37.73 

.19 

•9922u 

-.52 

.74 

1.00188 

40.65 

.20 

.99051 

-.59 

.75 

1.03995 

41.98 

.21 

.98857 

-.65 

.76 

1.08202 

43.21 

.22 

.98637 

-.72 

.77 

1.12839 

44.34 

.23 

.98389 

-.79 

.78 

1.17943 

45.36 

.24 

.98111 

-.86 

.79 

1.23555 

46.26 

.25 

.97802 

-.92 

.80 

1.29724 

47.03 

.26 

.97461 

-.99 

.81 

1.36506 

47.65 

.27 

.97087 

-1.04 

.82 

1.43964 

48.13 

.28 

.96679 

-1.09 

.83 

1.52174 

48.45 

.29 

.96235 

-1.12 

.84 

1.61222 

48.59 

.30 

.95756 

-1.15 

.85 

1.71205 

48.54 

.31 

.95241 

-1.16 

.86 

1.82236 

48.28 

.32 

.94691 

-1.16 

.87 

1.94438 

47.77 

.33 

.94106 

-1.13 

.88 

2.07945 

46.99 

.34 

.93486 

-1.09 

.89 

2.22898 

45.91 

.35 

.92832 

-1.02 

.90 

2.39429 

44.47 

.36 

.92147 

-.92 

.91 

2.57645 

42.63 

.37 

.91432 

-.79 

.92 

2.77585 

40.33 

.38 

.90688 

-.63 

.93 

2.99170 

37.49 

.39 

.89919 

-.44 

.94 

3.22103 

34.05 

.40 

.89127 

-.20 

.95 

3.45764 

29.94 

.41 

.88316 

•  08 

.96 

3.69080 

25.12 

.42 

.87490 

.40 

.97 

3.9C464 

19.60 

.43 

.86653 

.78 

.98 

4.07917 

13.45 

.44 

.85810 

1.21 

.99 

4.19390 

6.83 

.45 

.84966 

1.69 

1.00 

4.23378 

0.00 

.46 

.84127 

2.24 

1.01 

4.19468 

-6.77 

.47 

.83298 

2.85 

1.02 

4.08492 

-13.20 

.48 

.82487 

3.52 

1.03 

3.92177 

-19.09 

.49 

.81700 

4.27 

1.04 

3.72542 

-24.32 

.50 

.80946 

5.09 

1.05 

3.51397 

-28.86 

.51 

.80232 

5.98 

1.06 

3.30100 

-32.73 

.52 

.79566 

6.95 

1.07 

3.09533 

-35.99 

.53 

.78959 

7.99 

1.08 

2.90190 

-38.72 

.54 

.78420 

9.12 

1.09 

2.72300 

-40.97 

.55 

.77959 

10.32 

1.10 

2.55917 

-42.82 

.56 

.77586 

11.60 

1.11 

2.41000 

-44.33 

.57 

.77313 

12.95 

1.12 

2.27455 

-45.54 

.58 

.77150 

14.38 

1.13 

2.15166 

-46.50 

.59 

.77111 

15.87 

1.14 

2.04013 

-47.24 

.60 

.77206 

17.43 

1.15 

1.93881 

-47.80 

.61 

.77449 

19.04 

1.16 

1.84660 

-48.19 

.62 

.77852 

20. 70 

1.17 

1.76253 

-48.^5 

.63 

.78428 

22.40 

1.18 

1.68572 

-48.58 

.64 

.79191 

24.14 

1.19 

1.61542 

-48.60 
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C   C   ADPI (Z(S)-PHASE)  VS.  FREQUENCY,  OUTPUT  DATA,  TYPE  6  FILTER,  FIG. 42 


Wl 


2(  JW) 


1 


Wl 


Z(  JW) 


01 


.10 

.99976 

.12 

.99951 

.14 

.99912 

.15 

.99885 

.17 

.99814 

.19 

.99718 

.20 

.99659 

.22 

.99517 

.23 

.99434 

.24 

.99341 

.25 

.99240 

.26 

.99129 

.27 

.99009 

.28 

.98881 

.29 

.98743 

.30 

.98597 

.31 

.98443 

.32 

.98281 

.33 

.98112 

.34 

.97937 

.35 

.97757 

.36 

.97572 

.37 

.97384 

.38 

.97195 

.39 

.97004 

.40 

.96815 

.41 

.96629 

.42 

.96446 

.43 

.96269 

.44 

.96100 

.45 

•9594U 

.46 

.95792 

.47 

.95656 

.48 

.95536 

.49 

•95432 

.50 

.95347 

.51 

.95282 

.52 

.95239 

.53 

.95219 

.54 

.95223 

.55 

.95253 

.56 

.95308 

.57 

.95390 

.58 

.95498 

.59 

.95632 

.60 

.95791 

.61 

.95975 

.62 

.96180 

.63 

.96406 

.64 

.96648 

-.03 
-.05 
-.07 
-.09 
-.12 
-.15 
-.17 
-.20 
-.21 
-.22 
-.23 
-.24 
-.24 
-.24 
-.23 
-.22 
-.20 
-.18 
-.14 
-.10 
-.05 
.00 
.07 
.15 
.24 
.34 
.45 
.58 
.72 
.86 
1.02 
1.20 
1.38 
1.58 
1.78 
2.00 
2.23 
2.46 
2.70 
2.95 
3.20 
3.45 
3.71 
3.97 
4.22 
4.48 
4.72 
4.96 
5.20 
5.42 


.65 

.96903 

5.64 

.67 

.97437 

6.04 

.69 

.97964 

6.40 

.70 

.98210 

6.57 

.72 

.98634 

6.90 

.74 

.98922 

7.25 

.75 

.99001 

7.45 

.77 

.99013 

7.96 

.78 

.98949 

8.29 

.79 

.98846 

-351.15 

.80 

.98719 

-350.66 

.81 

.98588 

9.79 

.82 

.98486 

10.53 

.83 

.98456 

11.42 

.84 

.98559 

12.48 

.85 

.98871 

13.74 

.86 

.99494 

15.20 

.87 

1.00555 

16.86 

.88 

1.02216 

18.72 

.89 

1.04679 

20.73 

.90 

1.08194 

22.83 

.91 

1.13074 

24.93 

.92 

1.19704 

26.87 

.93 

1.28561 

28.47 

.94 

1.40212 

29.44 

.95 

1.55284 

29.46 

.96 

1.74272 

28.06 

.97 

1.97021 

24.70 

.98 

2.21421 

18.83 

.99 

2.41842 

10.32 

1.00 

2.50065 

0.00 

1.01 

2.41999 

-10.22 

1.02 

2.22346 

-18.55 

1.03 

1.99142 

-24.29 

1.04 

1.77552 

-27.69 

1.05 

1.59440 

-29.27 

1.06 

1.44918 

-29.58 

1.07 

1.33512 

-29.01 

1.08 

1.24649 

-27.87 

1.09 

1.17809 

-26.40 

1.10 

1.12564 

-24.74 

1.11 

1.08572 

-23.03 

1.12 

1.05563 

-21.33 

1.13 

1.03321 

-19.70 

1.14 

1.01679 

-18.18 

1.15 

1.00497 

-16.79 

1.16 

.99673 

-15.53 

1.17 

.991  18 

-14.40 

1.18 

.98767 

-13.40 

1.19 

.98566 

-12.53 
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C   ADPI (Z( S)-PHASE)  VS.  FREQUENCY,  OUTPUT  DATA,  TYPE  7  FILTER,  FIG. 45 


Wl 


Z(  JW) 


01 


Wl 


Z(  JW) 


01 


.10 

1  .00000 

.00 

.65 

.97784 

-2.36 

.12 

1.00000 

.00 

.67 

.96597 

-2.56 

.14 

1.00000 

.00 

.69 

.95132 

-2.65 

.15 

1.00000 

.00 

.70 

.94294 

-2.63 

.17 

1.00000 

.00 

.72 

.92414 

-2.43 

.19 

1.00002 

.00 

.74 

.90300 

-1.92 

.20 

1.00003 

.00 

.75 

.89179 

-1.53 

.22 

1.00008 

.00 

.77 

.86886 

-.33 

.23 

1.00011 

.01 

.78 

.85759 

.48 

.24 

1.00015 

.01 

.79 

.84682 

1.47 

.25 

1.00020 

.01 

.80 

.83693 

2.64 

.26 

1.00026 

.02 

.81 

.82836 

4.01 

.27 

1.00034 

.02 

.82 

.82165 

5.60 

.28 

1.00043 

.02 

.83 

.81744 

7.41 

.29 

1.00053 

.02 

.84 

.81648 

9.44 

.30 

1.000  66 

.02 

.85 

.81967 

11.68 

.31 

1.00080 

.02 

.86 

.82806 

14.09 

.32 

1.00096 

.02 

.87 

.84285 

16.64 

.33 

1.00115 

.02 

.88 

.86548 

19.27 

.34 

1.00135 

.02 

.89 

.89762 

21.88 

.35 

1.00158 

.01 

.90 

.94120 

24.36 

.36 

1.0C184 

.00 

.91 

.99854 

26.59 

.37 

1.00211 

-.00 

.92 

1.07233 

28.40 

.38 

1.00241 

-.01 

.93 

1.16570 

29.61 

.39 

1.00273 

-.03 

.94 

1.28203 

30.00 

.40 

1.00307 

-.04 

.95 

1.42418 

29.27 

.41 

1.00342 

-.06 

.96 

1.59257 

27. ]2 

.42 

1.00379 

-.08 

.97 

1.78068 

23.17 

.43 

1.00415 

-.11 

.98 

1.96766 

17.16 

.44 

1.00452 

-.15 

.99 

2.11307 

9.19 

.45 

1.00487 

-.18 

1.00 

2.16889 

90.08 

.46 

1.00521 

-.23 

1.01 

2.11415 

-9.11 

.47 

1.00551 

-.28 

1.02 

1.97447 

-16.88 

.48 

1.00576 

-.34 

1.03 

1.79753 

-22.72 

.49 

1.00595 

-.40 

1.04 

1.62057 

-26.63 

.50 

1.00606 

-.48 

1.05 

1.46201 

-28.90 

.51 

1.00607 

-.56 

1.06 

1.32734 

-29.89 

.53 

1.00568 

-.75 

1.08 

1.12512 

-29.20 

.54 

1 .00523 

-.85 

1.09 

1.05166 

-27.98 

.55 

1 .00457 

-.97 

1.10 

.99271 

-26.40 

.56 

1.00366 

-1.09 

1.11 

.94579 

-24.57 

.57 

1.00247 

-1.22 

1.12 

.90883 

-22.60 

.58 

1.00096 

-1.36 

1.13 

.88012 

-20.57 

.59 

.99908 

-1.50 

1.14 

.85826 

-18.53 

.60 

.99681 

-1.65 

1.15 

.84207 

-16.53 

.61 

.99408 

-1.80 

1.16 

.83056 

-14.61 

.62 

.99086 

-1  .95 

1.17 

.82291 

-12.79 

•  63 

.98711 

-2.09 

1.18 

.81842 

-11.09 

.64 

.98278 

-2.23 

1.19 

.81652 

-9.52 
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APPENDIX  F 

Voltage  Transfer  Function  (VTF)  Coefficients.   To  deter- 
mine the  VTF  frequency  response  of  the  previously  discussed 
filters,  it  is  necessary  to  know  the  VTF  coefficients. 
Since  all  of  the  filters,  with  the  exception  of  the  constont-1 
filter,  have  the  same  general  ladder  configuration  of  Fig.  56» 
a  pseudo-generalized  computer  program  was  written  to  deter- 
mine these  coefficients.   First  the  transfer  matrix  of  this 


f 


1/1 

■A/W 


*2 

AA/V 


*4 
■AA/V 


E. 
in 


1 


5         o 


1 


Fig.  56.   Generalized  ladder  configuration  of  filters 
being  analyzed. 


generalized  ladder  network  was  determined  and  the  inverse 
of  the  (1,1)  element  was  used  as  the  VTF, eC   .   The  general- 
ized expression  f  or  oC    is  given  in  eauation  (4-4).   Eauations 


cC   =  l+Y1+Y5+Y5+Z2Y5+Z2Y5+Z4Y5+Y1Z2Y$+yiZ4Y5 
+Y1^2Y5+Y3^4Y5+Z2Y3Z4Y5+Y1Z2Y3Z4Y5 


(44) 


Yn  = 


"nN 


and  Z   = 
n 


■SSL  8nd  ^    ^N 


■nD 


JnD 


cC 


*^*N  "  Y1DZ2DY3DZ4DY5D 


(45) 


92 


13 

n=l 


where 


Bl  =  Y1DZ2DY3DZ4DY5D    B7   =  Z4NY5NY1DZ2DY3D 

B2  =  Y1NZ2DY3DZ^DY5D    B8   =  Y1NZ2NY3NZ4DY5D 

B3  =  Y3NY1DZ2DZ4DY5D   B9   =  Y1NZ4NY5NZ2DY3D 

B4  =  Y5NY1DZ2DY3DZ4D   B10  =  Y1NZ2NY5NY3DZ4D     ^^ 
B5  =  Z2NY3NY1DZ4DY5D    Bll  =  Y3NZ4NY5NY1DZ2D 

B6  =  Z2NY5NY1DY3DZ4D    B12  =  Z2NY3NZ4NY5NY1D 

B13  =  Y1NZ2NY3NZ4NY5N 

(45),  (46),  and  (47)  put  j£    in  the  form  which  the  computer 
program  can  handle  it.   The  following  restrictions  were 
placed  on  ladder  component  expressions  to  simplify  the  com- 
puter program: 

1.  Zpr)  =  Yxn  =  Z^c  =  1  for  all  types. 

2.  ZpN  Y^«,  and  Z^  are  of  the  form  -  As. 

3.  Y-jN  and  Y,~  are  a  maximum  of  fifth  order  polynomials. 

4.  YcN  and  Y,-*   are  a  maximum  of  second  order  polynom- 
ials. 

Using  these  restrictions,  oC~   reduces  to  equation  (48)  and 
cCn   reduces  to  equation  (49).   With  these  restrictions, 
Type  5  exceeds  the  restrictions  on  Y,^  and  Y,~  so  it  was 
ignored. 


95 


-c  - 


D 


Y1DY5D(1+Y5N+Z2Y3N)   +   Y1NY5D(1+Z2NY3N> 


+   Y 


1DY5N^1+Z2N+Z4N+Y3NZ4N",'Z2NY3NZ4N^  (48) 


+    Y1NY5N(Z2N+Z4N+Z2NY3NZ4N) 


^"N   =    Y1DY5D 


(49) 
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APPENDIX  G 

FORGO  Programs  for  VTF  Coefficient  Calculations.   The 
computer  program  on  page?  94  through  97  was  written  to 
evaluate  the  coefficients  of  eouations  (45),  (48)  and  (49). 
The  input  and  output  data  for  all  filters  except  the  constant-1 
filter  and  Type  5  are  given  on  the  pages  following  the  com- 
puter program. 


READ  Y1N(I),  Y1D(I) 
READ  Y5N(I),  Y5D(I) 
READ  Zow,  Y™,  Z. 


'2N 


3N'   4N 


I 


Y1DY5D=E1D5D(I) 
Y1NY5D-B1N5D(I) 

Y1NY5N=E1N5N(I) 
Y1DY5N=E1D5N(I) 


SL 


A(1)-E1D5D(I) 


13 


B(I)  =  ]T  Bn(l) 


n=l 


I 


i 


PUNCH 
COEFFIC. 
IDENTIF- 
ICA. 


JSZ. 


PUNCH 
A(I), 

B(I) 


Fig.  57.   Block  diagram  of  VTF  coefficient  computer  progrsm. 


c  c 


V 

D 
D 
D 
D 
18  R 
R 
R 
E 
F 
F 
F 
E 
E 
E 
E 
E 
E 
F 
F 
F 
F 
F 
E 
E 
E 
E 
F 
E 
E 
F 
F 
F 
F 
F 
E 
E 
F 
E 
E 
F 
F 
G 
P 
F 
F 
A 
D 

1  A 
A 
A 
A 
B 
D 

2  B 


OLTAGt  TRANSFE 

IMENSICN  E1N5D 

IMENSICN  B2( 10 

IMENSICN  B10(l 

IMENSICN  Y5N(3 

FAD,Y1N0, ( Y1N( 

FAD,Z2N,Y3N,Z4 

FAD,Y5N0,Y5N( 1 

1ND0  =Y1N0*Y5D 

1N5FM 1 )=Y1N( 1 ) 

1N5D(2)=Y1N(2) 

1N5D(3)=Y1N(3) 

1N5D(4)=Y1N(4) 

1N5D(5)=Y1N( 5) 

1N5D(6  )  = 

1N5D(7)= 

1NN0  =Y1N0*Y5N 

1N5NI 1 )=Y1N( 1 ) 

1N5N(2)=Y1N(2) 

1N5N(3 )=Y1N(3 ) 

1N5N(4)=Y1N(4) 

1N5N(5)=Y1N( 5) 

1N5N(6 )= 

1N5N(7)= 

1DN0  =Y1DC*Y5N 

1D5N(1 )=Y1D( 1 ) 

1D5N(2)=Y1D(2) 

1D5N(3)=Y1D<3) 

1D5N(4)=Y1D(4) 

1D5N(5)=Y1D( 5) 

1D5N(6)= 

1D5N(7)= 

1DD0  =  > 1D0*Y5D 

1D5D(1 )=Y1D(1 ) 

1D5D(2)=Y1D(2) 

1D5D(3)=Y1D(3) 

1D5D(4)=Y1D(4) 

1D5D(5)=Y1D(5) 

1D5D(6)= 

1D5D(7)= 

=  Y3N 

1=Z2N*Y3N 

=  Z2N 

=  Z4N 

2=Y3N*Z4N 

3=F2*G 

10=E1DOO 

C  1  1=1*7 

1 (  I  )=E1D5D(  I  ) 

1  (  8  )  =  0  •  0 

1 (9)=0.0 
1 ( 10)=0.0 

102asA10 

C  2  1  =  1, 1C 
1  (  I  )=A1(  I  ) 


ICN(VTF)  CCE 

N5N( 10 ) ,E1D5 

) »B4( 10) »B5( 

10) ,B12 (10)  , 

) 

5)  ,Y1D0»  (Y1D 


( I ) ,1=1,5) 
)  ,Y5D0,Y5D( 1 ) ,Y5D(2  ) 


1N0*Y5D( 1  ) 
1N( 1 )*Y5D(  1  ) 
1N(2 )*Y5D( 1  ) 
1N(3)*Y5D( 1 ) 
1N(4)*Y5D(  1  ) 
1N(5 )*Y5D( 1 ) 


R  FUNCT 
( 10) ,E1 
),B3(  10 
0) ,B11  ( 
) ,YSD(3 
I  )  ,  1  =  1  , 
N 

) ,Y5N  (2 
0 

♦Y5D0+Y 
*Y5D0+Y 
*Y5D0+Y 
*Y5D0+Y 
*Y5D0+Y 
Y 

0 

*Y5NO+Y1NO*Y5N( 1  ) 
*Y5N0+Y1N( 1 )*Y5N(  1  ) 
*Y5N0+Y1N(2 )*Y5N( 1  ) 
*Y5N0+Y1N(3 )*Y5N( 1  ) 
*Y5N0+Y1N(4)*Y5N( 1 ) 
Y1N(5 )*Y5N<1  ) 

0 

*Y5N0+Y1D0*Y5N( 1  ) 
*Y5N0+Y1D( 1 )*Y5N(1  ) 
*Y5N0+Y1D(2 )*Y5N(1  ) 
*Y5N0+Y1D(3 )*Y5N( 1  ) 
*Y5N0+Y1D(4)*Y5N(  1  ) 
Y1D<5)*Y5N( 1  ) 

0 

*Y5D0+Y1D0*Y5D( ] ) 
*Y5D0+Y1D( 1 )*Y5D( 1 ) 
*Y5D0+Y1D(2)*Y5D<1) 
*Y5D0+Y1D(3 )*Y5D(1 ) 
*Y5D0+Y1D(4)*Y5D( 1 ) 
Y1D(5)*Y5D(1) 


95 

FFICIENT  EVALUATION 
N(10),E1D5D(10),A1(10),B1(10) 
10),B6(10),B7(10),B8(10),B9(10 
D13(10),B(10),Y1N(5),Y1D(5) 


+Y1N0*Y5D(2) 
+Y1N(1 )*Y5D( 2) 
+Y1N(2)*Y5D( 2) 
+Y1N(3 )*Y5D( 2) 
+Y1N(4)*Y5D( 2) 
Y1N(5)*Y5D( 2) 


+Y1N0*Y5N(2) 
+Y1N( 1 )*Y5N( 2 ) 
+Y1N(2)*Y5N(2) 
+Y1N(3)*Y5N(2) 
+Y1N(4)*Y5N(2) 
Y1N( 5)*Y5N(2) 


+Y1D0*Y5N(2) 
+  Y1D( 1 )*Y5N(2) 
+Y1D(2)*Y5N(2) 
+Y1D(3)*Y5N(2) 
+Y1D(4)*Y5N(2) 
Y1D(5)*Y5N(2) 


+Y1D0*Y5D(2) 
+  Y1D( 1)*Y5D(2) 
+Y1D(2)*Y5D(2 ) 
+Y1D(3)*Y5D( 2 ) 
+Y1D(4)*Y5D(2) 
Y1D(5)*Y5D(2) 
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11 


B20=E1N00 

DC 

3     1=1,7 

B2I 

I )=E1N5D( I ) 

B2I 

8  )  =  0  .  0 

B2I 

9)=0.0 

B2I 

10)=0.0 

B3C=0.U 

B3  ( 

1 )=F*E1DD0 

DC 

4    1=2,8 

B3I 

I )=F*E1D5D( 1-1 ) 

B3< 

9)=0.0 

B3( 

10)=0.0 

B4C=E1DN0 

DC 

5    1=1,7 

B4< 

I  )=E1D5N(  I  ) 

B4( 

8)=o.O 

B4< 

9)=0.0 

B4( 

10)=0.0 

650=0. u 

B5( 

1)=0.0 

B5I 

2)=F1*E1DDG 

DC 

6    1=3,9 

B5< 

I )=F1*E1D5D( 1-2) 

B5( 

10)=0.0 

B6C 

,  =  0.0 

B6( 

1 )=G*E1DN0 

DC 

7     1=2,8 

B6( 

I )=G*E1D5N( 1-1 ) 

B6( 

9)=0.0 

B6( 

10)=0.0 

B70=G.O 

B7( 

1 )=P*E1DN0 

DC 

8    1=2,8 

B7< 

I )=P*E1D5N< 1-1 ) 

B7( 

9)=0.0 

B7( 

10)=0.0 

B80=0.0 

B8( 

1 )  =  0  .  o 

B8I 

2 )=F1*E1ND0 

DC 

9     1=3,9 

B8( 

I )=F]*E1N5D( 1-2) 

B8< 

10)=0.0 

B90=0.c 

B9I 

1 )=P*E1NN0 

DC 

10    1=2,8 

B9( 

I )=P*E1N5N( 1-1 ) 

B9( 

9)=u.O 

B9I 

10)=0.0 

B100=0.0 

B10(1)=G*E1NN0 

DC 

11     1=2,8 

B10( I )=G*E1N5N( 1-1 ) 

BlCi  (9)  =0.0 

B10( 10)=0.C 

B110=0.0 

97 


12 


13 


Bll  (  1  )=0.0 
Bll  (2)=F2*E1DN0 
DO    12     1=3,9 
Bll  (  I  )=F2*E1D5N( 1-2) 
Bll (  10)=  J.O 
6120=0.0 
B12(  1  )=0.0 
B12(2)=0.0 
B12(3)=F3*E1DN0 
DC    13    1=4, lu 
B12(  I  )=F3*E1D5N<  1-3) 
B130=0.0 
B13(  1)=0.0 
B13(2)=0.0 
B13(3)=F3*E1NN0 
DC    14    1=4,10 
B13(  I  )=F3*E1N5N( 1-3) 
B0=B10Z+B20+B40 
DC    19    1=1,10 

BT    =    Bl ( I  )+B2(  I  )+B3(  I  )+B4(  I  )+B5 (  I  )+B6(  I  ) +B7 (  I  )+B8(  I  )+B9(  I  )+BlO( I 
19    B(I)=BT    +B1K  I  )+B12(  I  )+B13(  I  ) 
PUNCH    15,A10,B0 

FCRMAT(5X,10H    A(     0)        =    F10.5,5X,10H    B(     0) 
DC    16    1=1 ,10 
PUNCH    17,  I ,A1(  I  ) ,1 ,B( I  ) 

FCRMAT(5X,3H    A(I2»5H)        =  F10.5,5X,3H    B(I2»5H) 

GC    TC    18 
END 


14 


15 

16 
17 


=    F10.5) 


F10.5) 
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VTF  COEFFICIENT  EVALUATION,  INPUT  DATA,  ZOBEL  FILTER,  M=0.6 

.0   .6   .0   .0   .0   .0     1.   .0   .64   .0   .0   .0 

1.6   2.   1.6 

1.   .6   .64     1.   .0   .64 


VTF  COEFFICIENT  EVALUATION,  INPUT  DATA,  ZOBEL  FILTER,  M=0.707 
.0   .70710678   .0   .0   .0   .0     1.   .0   .5   .0   .0   .0 
1.70710678   2.   1.70710678 
1.   .70710678   .5     1.   .0   .5 


VTF  COEFFICIENT  EVALUATION,  INPUT  DATA,  TYPE  1  FILTER 

.0  2.   2.   1.   .0   .0     2.   2.   2.   1.   .0   .0 

1.  2.   1. 

1.  .0   .0     1.   .0   .0 


VTF  COEFFICIENT  EVALUATION,  INPUT  DATA,  TYPE  2  FILTER 

.o       2.       2.       1.       .0       .0  3.236       3.236       2.618       1.       .0       .0 

1.618       2.       1. 

1.       .0       .0  1.       .0       .0 


VTF  COEFFICIENT  EVALUATION,  INPUT  DATA,  TYPE  3  FILTER 
.0  2.155  2.155  1.  .0  .0  3.7325  3.7325  2.732 
1.732   2.155   1. 

1  .   .0   .0     1.   .0   .c 


1. 


VTF  COEFFICIENT  EVALUATION,  INPUT  DATA,  TYPE  4  FILTER 

•o   .665   1.33   1.   .0   .0     .9975   1.995   2.   1.   .0   .0 

1.5   1.333   .5 

1  .   .0   .0     1 .   .0   .0 


VTF  COEFFICIENT  EVALUATION,  INPUT  DATA,  TYPE  6  FILTER 
.0   2.0         3.2360679   3.85^+1019   2.61803399   1. 
3.2360679   5.2360679   7.23606789  5.8541019   3.2360679 
1.61803399   2.   1.61803399 
1.   .61803399   .61803399     1.   .0   .61803399 


VTF  COEFFICIENT  EVALUATION,  INPUT  DATA,  TYPE  7  FILTER 

.0   .8395616   1.7143309   1.9662335   1.3674296   .51436894 

1.3247343   2.7050226   3.5343409   3.0394495   1.6934152   .51436894 

1.5778883   1.9441298   1.7143309 

1.   .7143309   .4897314     1.   .0   .4897314 
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VTF  CCtFFICIENT  EVALUATION,  OUTPUT  DATA,  ZOBEL  FILTER,  FIG. 22 
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A(  1  ) 

A(  2) 

A(  3) 
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A(  9) 
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0.00000 
1.28000 
0.00000 
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0.00000 
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B( 

3) 

B( 

4) 
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5) 

B( 

6) 

B( 

7) 

B( 

8) 

B( 

9) 

B(  10) 


2.00000 

6.40000 

12.80000 

17.53600 

17.61280 

12.98432 

6.55360 

2.09715 

0.00000 

0.00000 

0.00000 


C   VTF  COEFFICIENT  EVALUATION,  OUTPUT  DATA,  ZOBEL  FILTER,  FIG. 25 
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1.00000 
0.00000 
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0.00000 
0.00000 
0.00000 
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B( 
B( 
B( 
B( 
B( 
B( 
B( 
B( 
B( 


0) 
]  ) 
2) 
3) 
4) 
5) 
6) 
7) 
8) 
9) 


B(  10) 


2.00000 

6.82843 

13.65685 

18.48528 

17.98528 

12.51041 

5.82843 

1.45711 

0.00000 

0.00000 

0.00000 


C   VTF  COEFrICIENT  EVALUATION,  OUTPUT  DATA,  TYPE  1  FILTER,  FIG. 28 
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C   C   VTF  COEFFICIENT  EVALUATION*  OUTPUT  DATA,  TYPE  2  FILTER,  FIG. 31 
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A{ 
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At 

7) 
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C   C   VTF  COEFFICIENT  EVALUATION,  OUTPUT  DATA,  TYPE  3  FILTER,  FIG. 34 
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C   C   VTF  COEFFICIENT  EVALUATION,  OUTPUT  DATA,  TYPE  4  FILTER,  FIG. 37 
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C   VTF  COEFFICIENT  EVALUATION,  OUTPUT  DATA,  TYPE  6  FILTER,  FIG. 42 
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C   C   VTF  COEFFICIENT  EVALUATION,  OUTPUT  DATA,  TYPE  7  FILTER,  FIG. 45 
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40.53567 

79.52404 

116.33964 
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APPENDIX  H 

FORGO  Programs  for  Sinusoidal  Response  Calculations 
of  the  VTF.   To  make  possible  the  computation  of  the  sin- 
usoidal VTF  response,  the  FORGO  digital  computer  program 
shown  on  pages  102  and  103  was  written  and  used.   The  input 
and  output  data  for  each  of  the  previously  discussed  filter 
configurations,  except  Type  5*  are  given  on  the  pages  fol- 
lowing the  computer  program. 


102 


READ  A(I) 
READ  B(I) 


I 


K=0 


( oG:=io,]20 


I 


K=K+1 

H=I 

V/=H/100. 

n(K)=w 


I 


T(K)  = 


T(s) 


;=0W1(K) 


T 


THE(K)=ARG  T(K) 


PUNCH 

COLUMN 
HWA  DINGS 


i 


K=l,110 


I 


PUNCH 
W1(K) 
T(K) 

THE(K) 


Fig.  58»   Block  diagram  of  VTF  computer  program. 
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4 

5 

45 


8 
11 
13 

12 

10 
16 

15 

14 

9 

28 
21 
23 

22 

19 
25 

24 

20 
2 


VTF( 

DIMEN 

FORMA 

FORMA 

READ, 

K  =  0 

DO  2 

K  =  K  +  1 

H=I 

W  =  H/1 

WKK) 

C  =  AO- 

D  =  A(  1 

E  =  BO- 

F  =  P(1 

T(K)  = 

IF(AB 

THF1  = 

GO  TO 

IF(C) 

IF(D) 

THE1  = 

GO  TO 

THE1 

GO  T 

!F(D 

THE1 

GO  T 

THE1 

IF(A 

THE2 

GO  T 

IF(E 

IF(F 

THE2 

GO  T 

THE2 

GO  T 

IF(F 

THE2 

GO  T 

THE2 

THE( 

CONT 

PUNC 

DO  7 

PUNC 

CONT 

GO  T 

FND 


T(S)  -  PHASE  SHIFT  )  VS.  FREQUENCY,  ALL  FILTERS 

SI  ON  A( 7) ,B( 10)  •W1(112)»T(112)»THE(112) 

T( /7X,2HW1 ,5X,5HT( JW) ,6X»2H01,9X,2HW1,5X,5HT( JW) ♦6X,2H01//) 

T(  2(5X,F5.2,3X,F7.5,2X,F7.2)  ) 

A0,(A(I),I=1,7),B0,(B(I)»I=1,10) 

1=10,120 


00, 

=  W 

A( 2) *W* 

)*W-A( 3 

B(2)*W* 

)*W-B( 3 

SQRT( (C 

S(C)-.l 

0. 

14 
10»10»1 
12»12.1 
ATAN(D/ 

14 


*2+A(4)*W**4-A (6  )*W**6 

)*W**3+A ( 5)*W**5-A( 7)*W**7 

*2+B(4)*W**4-B(6 )*W**6+B(8 )*W**8-B( 10 )*W**10 

)*W**3+B( 5)*W**5-B(7)*W**7+B( 9)*W**9 

** 2+0**2 )/ (E**2+F**2 ) ) 

F-40)3,3,8 


0*57.245779 
0*57.24779 


=ATAN (D/i 

0  14 

)15»15.16 

=180. -ATAN(D/ABS(C))*5 7. 245779 

0  14 

=  180.  +  AT, 

BS(E)-.l 

=  0. 

0  20 

)19,19,2 

)22»22»2 

=ATAN(F/I 

0  20 

=  ATAN(F/I 

0  20 

)24*24,25 

=180.-ATAN(F/APS(E))*57.245  779 

0  20 

=l80.+ATi 

K)=THE1-' 

INUE 

H  4 

J=l »55 
H  5»Vl(J),T{J)>THE(J),Wl(J+55)»T(J+55)»THE(J+55) 
INUE 
0  45 


AN(D/C)*57. 245779 
E-40)9,9,28 


1 
3 
E)*57. 245779 

E)*57. 245779 


AN(F/E)*57. 245779 
THF2 
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VTF  COEFFICIENTS,  INPUT  DATA,  CCNSTANT-1  FILTER 

1.   .0   .0   ,0   .0   ,0   .0   .0 

2  .   4.   4  .   2.   • 0   .0   .0   .0   • 0   .0   • C 


VTF  COEFFICIENTS,  INPUT  3ATA,  ZOBEL  FILTER,  M=0.6 

1.  .0   1.28   .0   .4096   .0   .0   .0 

2.  6.4   12.8   17.536   17.6128   12.98432   6.5536   2.09715   .0   .0   .0 

VTF  COEFFICIENTS,  INPUT  DATA,  ZOBEL  FILTER,  M=0.707 

1.  .0   1.   .0   .25   .0   .0   .0 

2.  6.82843   13.65685   18.48528   17.98528   12.51041   5.82843   1.45711 

,n   .0   .0 


VTF  COEFFICIENTS,  INPUT  DATA,  TYPE  1  FILTER 

2.   2.   2.   1  .   .0   .0   .0   .0 

4.   14.   26.   31.   26.   14.   4.   .0   .0   .0   .0 


VTF  COEFFICIENTS,  INPUT  DATA,  TYPE  2  FILTER 

3.236   3.236   2.618   1.   .0   .0   .0  .0 

6.472   23.41585   44.35954   54.21332  44.35954   23.41585   6.472 

•  o   .0 


0   .0 


VTF  COEFFICIENTS,  INPUT  DATA,  TYPE  3  FILTER 

3.7325   3.7325   2.732   1.   .0   .0   .0   .0 

7.465   27.86^73   53.72213   66.18855   53.72185   27.86045   7.46492 

.0   .0   •  0 


.0 


VTF  COEFFICIENTS,  INPUT  DATA,  TYPE  4  FILTER 
.9975   1.995   2.   1.   .0   .0   .0   .0 
1.995   7.97967   15.96867   19.97159   15.98367 
.0   .0 


7.99467   1.9995   .0   .0 


VTF  COEFFICIENTS,  INPUT  DATA,  TYPE  6  FILTER 

3.2360679  5.2360679  9.2360679  9.0901699  7.7082039  4.6180339  2.  .61803398 
6.4721359  31.416407  86.249223  164.82468  237.74767  267.83784  237.74767 
164.82468  86.249223  31.416407  6.4721359 


VTF  COEFFICIENTS,  INPUT  DATA,  TYPE  7  FILTER 

1.85451084       3.78679208      5.85597666 

4.79370693 

2.80385847        1.16097257      0.35264130 
3.7090216        19.7844959      56.7463506 

162.8651935 

184.2725252      162.8651935     111.3266175 

19. 7844959 

3.7090216 


6.10947126 


111.3266175 


56.7463506 
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C   C   VTF(T(S)-PHASE)  VS.  FREQUENCY*  OUTPUT  DATA,  CONSTANT-1  FILTER,  FIG. 19 
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.48054 
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-84.31 

.14 
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.49999 
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.49996 

-26.58 

.78 

.45172 

-101.36 

.24 

.49995 

-27.76 

.79 

.44846 

-102.93 

.25 

.49994 

-28.94 

.80 

.44506 

-104.50 

.26 

.49992 

-30.13 

.81 

.44152 

-106.07 

.27 

.49990 

-31.32 

.82 

.43785 

-107.64 

.28 

.49988 

-32.51 

.83 

.43405 

-109.21 

.29 

.49985 

-33.71 

.84 

.43012 

-110.78 

.30 

.49982 

-34.91 

.85 

.42607 

-112.34 

.31 

.49978 

-36.12 

.86 

.42189 

-113.91 

.32 

.49973 

-37.33 

.87 

.41759 

-115.47 

.33 

.49968 

-38.55 

.88 

.41318 

-117.02 

.34 

.49961 

-39.77 

.89 

.40866 

-118.57 

.35 

.49954 

-41.00 

.90 

.40404 

-120.12 

.36 

.49946 

-42.23 

.91 

.39931 

-121.66 

.37 

.49936 

-43.47 

.92 

.39450 

-123.19 

.38 

.49925 

-44.72 

.93 

.38960 

-124.71 

.39 

.49912 

-45.97 

.94 

.38463 

-126.22 

,40 

.49898 

-47.22 

.95 

.37958 

-127.72 

.41 

.49882 

-'+  8.  49 

.96 

.37448 

-129.21 

.42 

.49863 

-+9.76 

.97 

.36931 

-130.68 

.43 

.49843 

-51.04 

.98 

.36410 

-132.15 

.44 

.49820 

-52.32 

.99 

.35884 

-133.60 

.45 

.49794 

-53.62 

1.00 

.35355 

-135.04 

.46 

.49765 

-54.92 

1.01 

.34824 

-136.46 

.47 

.49733 

-56.23 

1.02 

.34290 

-137.87 

.48 

.49697 

-57.54 

1.03 

.33756 

-139.27 

.49 

.49658 

-58.87 

1.04 

.33220 

-140.64 

.50 

.49614 

-60.20 

1.05 

.32685 

-142.01 

.51 

.49566 

-61.55 

1.06 

.32151 

-143.35 

.52 

.49513 

-62.90 

1.07 

.31618 

-144.68 

.53 

.49455 

-64.26 

1.08 

.31087 

-145.99 

.54 

.49391 

-65.63 

1.09 

.30559 

-147.29 

.55 

.49322 

-67.01 

1.1C 

.30034 

-148.56 

.56 

.49246 

-68.40 

1.11 

.29512 

-149.82 

.57 

.49164 

-69.80 

1.12 

.28994 

-151.06 

.58 

.49075 

-71.21 

1.13 

.28481 

-152.29 

.59 

.48978 

-72.62 

1.14 

.27973 

-153.49 

.60 

.48873 

-74.05 

1.15 

.27470 

-154.68 

.61 

.48760 

-75.49 

1.16 

.26972 

-155.85 

.62 

.48638 

-76.94 

1.17 

.26481 

-157.00 

.63 

.48507 

-78.39 

1.18 

.25995 

-158.13 

.64 

.48366 

-79.86 

1.19 

.25516 

-159.25 
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C   C   VTF(T(S)-PHASE)  VS.  FREQUENCY*  OUTPUT  DATA,  ZOBEL  FILTER,  FIG. 22 


Wl 


T(JW) 


01 


Wl 


T(  JW) 


01 


.10 

.50000 

-18.36 

.65 

.49980 

-135.43 

.11 

.50000 

-20.21 

.66 

.49981 

-138.20 

.12 

.50000 

-22.06 

.67 

.49984 

-141.02 

.14 

•  5000o 

-25.77 

.69 

.49988 

-146.80 

.15 

.50000 

-27.63 

.70 

.49990 

-149.78 

.17 

.50000 

-31.37 

.72 

.49994 

-155.92 

.19 

.50000 

-35.12 

.74 

.49997 

-162.33 

.20 

•  50000 

-37.01 

.75 

.49998 

-165.64 

.22 

.50000 

-40.80 

.77 

.50000 

-172.53 

.23 

.49999 

-42.70 

.78 

.50000 

-176.11 

.24 

.49999 

-44.61 

.79 

.50000 

-179.79 

.25 

.49999 

-46.53 

.80 

.50000 

-183.58 

.26 

.49999 

-48.45 

.81 

.50000 

-187.49 

.27 

.49999 

-50.39 

.82 

.50000 

-191.52 

.28 

.49998 

-52.33 

.83 

.50000 

-195.70 

.29 

.49998 

-54.28 

.84 

.50000 

-200.04 

.30 

.49997 

-56.24 

.85 

.49998 

-204.55 

.31 

.49997 

-58.21 

.86 

.49994 

-209.26 

.32 

.49996 

-60.19 

.87 

.49984 

-214.18 

.33 

.49996 

-62.18 

.88 

.49964 

-219.35 

.34 

.49995 

-64.18 

.89 

.49924 

-224.81 

.35 

.49994 

-56.19 

.90 

.49852 

-230.59 

.36 

.49994 

-68.21 

.91 

.49727 

-236.74 

.37 

.49993 

-70.25 

.92 

.49516 

-243.33 

.38 

.49992 

-72.30 

.93 

.49174 

-250.40 

.39 

.49991 

-74.36 

.94 

.48635 

-258.03 

.40 

.49990 

-76.44 

.95 

.47811 

-266.27 

.41 

.49989 

-78.53 

.96 

.46595 

84.68 

.42 

.49987 

-80.64 

.97 

.44873 

75.15 

.43 

.49986 

-82.76 

.98 

.42551 

65.05 

.44 

.49985 

-84.90 

.99 

.39596 

54.52 

.45 

.49984 

-87.06 

1.00 

.36070 

43.79 

.46 

.49982 

-89.23 

1.01 

.32136 

33.14 

.47 

.49981 

-91.58 

1.02 

.28030 

22.81 

.48 

.49980 

-93.79 

1.03 

.23990 

13.03 

.49 

.49978 

-96.03 

1.04 

.20210 

3.90 

.50 

.49977 

-98.28 

1.05 

.16809 

-4.52 

.51 

.49976 

-100.56 

1.06 

.13837 

-12.24 

.52 

.49975 

-102.86 

1.07 

.11293 

-19.30 

.53 

.49974 

-105.18 

1.08 

.09149 

-25.78 

.54 

.49973 

-107.53 

1.09 

.07361 

-31.74 

.55 

.49973 

-109.91 

1.10 

.05880 

-37.24 

.56 

.49972 

-112.31 

1.11 

.04661 

-42.34 

.57 

.49972 

-114.74 

1.12 

.03663 

-47.09 

.58 

.49972 

-117.20 

1.13 

.02850 

-51.53 

.59 

.49973 

-119.70 

1.14 

.02191 

-55.70 

.60 

•4°973 

-122.23 

1.15 

.01659 

-59.62 

.61 

.49974 

-124.79 

1*16 

.01234 

-63.34 

.62 

.49975 

-127.39 

1.17 

.00896 

-66.86 

.63 

.49976 

-130.03 

1.18 

.00632 

-70.21 

.64 

.49978 

-132.71 

1.19 

.00429 

-73.40 
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C   C   VTFCT(S)-PHASE)  VS.  FREQUENCY,  OUTPUT  DATA,  ZOBEL  FILTFR,  FIG. 25 


Wl 


T(  JW) 


Wl 


T(JW) 


01 


.10 

.50000 

-19.59 

.65 

.49988 

-143.43 

.11 

.50000 

-21.56 

.66 

.49987 

-146.31 

.12 

.50000 

-23.53 

.67 

.49987 

-149.23 

.14 

.50000 

-27.49 

.69 

.49988 

-155.22 

.15 

.50000 

-29.47 

.70 

.49989 

-158.29 

.17 

.50000 

-33.46 

.72 

.49993 

-164.62 

.19 

.50000 

-37.46 

.74 

.49997 

-171.21 

.20 

.50000 

-39.47 

.75 

.49999 

-174.61 

.22 

.50000 

-43.50 

.77 

.50000 

-181.67 

.23 

.50000 

-45.53 

.78 

.49998 

-185.33 

.24 

.50000 

-47.56 

.79 

.49992 

-189.10 

.25 

•  50000 

-49.61 

.80 

.49981 

-192.97 

.26 

.50000 

-51.65 

.81 

.49961 

-196.97 

.27 

•  50000 

-53.71 

.82 

.49930 

-201.10 

.28 

.50000 

-55.78 

.83 

.49882 

-205.38 

•  29 

.50000 

-57.85 

.84 

.49811 

-209.81 

•  30 

.50000 

-59.94 

.85 

.49709 

-214.42 

.31 

.50000 

-62.03 

.86 

.49565 

-219.21 

.32 

.50000 

-64.13 

.87 

.49366 

-224.22 

.33 

.50000 

-66.25 

.88 

.49093 

-229.44 

.34 

.50000 

-68.37 

.89 

.48727 

-234.91 

.35 

.50000 

-70.51 

.90 

.48242 

-240.63 

.36 

.50000 

-72.66 

.91 

.47609 

-246.61 

.37 

.50000 

-74.82 

.92 

.46799 

-252.86 

.38 

.50000 

-77.00 

.93 

.45780 

-259.38 

.39 

.50000 

-79.18 

.94 

.44523 

-266.15 

.40 

.50000 

-81.39 

.95 

.43009 

86.70 

.41 

.50000 

-83.60 

.96 

.41228 

79.53 

.42 

.49999 

-85.84 

.97 

.39191 

72.23 

.43 

.49999 

-88.08 

.98 

.36924 

64.88 

.44 

.49999 

-90.50 

.99 

.34473 

57.55 

.45 

.49999 

-92.79 

1.00 

.31898 

50.32 

.46 

.49999 

-95.09 

1.01 

.29264 

43.24 

.47 

.49998 

-97.41 

1.02 

.26638 

36.38 

.48 

.49998 

-99.74 

1.03 

.24077 

29.78 

.49 

.49998 

-102.10 

1.04 

.21629 

23.46 

.50 

.49997 

-104.48 

1.05 

.19328 

17.45 

.51 

.49997 

-106.88 

1.06 

.17194 

11.73 

.52 

.49997 

-109.31 

1.07 

.15237 

6.32 

.53 

.49996 

-111.75 

1.08 

.13459 

1.19 

.54 

.49995 

-114.23 

1.09 

.11856 

-3.68 

.55 

.49995 

-116.73 

1.10 

.10418 

-8.29 

.56 

.49994 

-119.25 

1.11 

.09133 

-12.66 

.57 

.49993 

-121.81 

1.12 

.07991 

-16.82 

.58 

.49993 

-124.39 

1.13 

.06976 

-20.78 

.59 

.49992 

-127.00 

1.14 

.06079 

-24.56 

.60 

.49991 

-129.65 

1.15 

.05285 

-28.17 

.61 

.49990 

-132.33 

1.16 

.04585 

-31.62 

.62 

.49989 

-135.05 

1.17 

.03968 

-34.94 

.63 

.49989 

-137.80 

1.18 

.03424 

-38.12 

.64 

.49988 

-140.59 

1.19 

.02947 

-41.17 
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C   VTF(T(S)-PHASE)  VS.  FREQUENCY,  OUTPUT  DATA,  TYPE  1  FILTER,  FIG. 28 


Wl 


T(  JW) 


01 


Wl 


T(  JW) 


01 


.10 

.49935 

-14.34 

.65 

.44796 

-98.59 

.11 

.49920 

-15.78 

.66 

.44641 

-100.23 

.12 

.49904 

-17.23 

.67 

.44488 

-101.89 

.14 

.49866 

-20.12 

.69 

.44192 

-105.23 

.15 

.49845 

-21.57 

.70 

.44051 

-106.93 

.17 

.49795 

-24.48 

.72 

.43781 

-110.36 

.19 

.49737 

-27.40 

.74 

.43532 

-113.87 

.20 

.49704 

-28.86 

.75 

.43417 

-115.66 

.22 

.49631 

-31.80 

.77 

.43202 

-119.34 

.23 

.49590 

-33.27 

.78 

.43102 

-121.23 

.24 

.49547 

-34.75 

.79 

.43007 

-123.16 

.25 

.49501 

-36.23 

.80 

.42914 

-125.14 

.26 

.49451 

-37.71 

.81 

.42824 

-127.16 

.27 

.49398 

-39.20 

.82 

.42732 

-129.25 

.28 

.49342 

-40.69 

.83 

.42638 

-131.39 

.29 

.49283 

-42.19 

.84 

.42536 

-133.61 

.30 

.49220 

-43.68 

.85 

.42422 

-135.90 

.31 

.49154 

-45.19 

.86 

.42290 

-138.28 

.32 

.49084 

-46.69 

.87 

.42134 

-140.74 

.33 

.49010 

-48.20 

.88 

.41945 

-143.30 

.34 

.48932 

-49.72 

.89 

.41713 

-145.96 

.35 

.48850 

-51.23 

.90 

.41427 

-148.72 

.36 

.48764 

-52.75 

.91 

.41075 

-151.58 

.37 

.48675 

-54.28 

.92 

.40644 

-154.54 

.38 

.48581 

-55.81 

.93 

.40121 

-157.60 

.39 

.48483 

-57.34 

.94 

.39496 

-160.74 

.40 

.48381 

-58.88 

.95 

.38758 

-163.94 

.41 

.48275 

-60.42 

.96 

.37900 

-167.19 

.42 

.48164 

-61.96 

.97 

.36923 

-170.45 

.43 

.48050 

-63.50 

.98 

.35828 

-173.70 

.44 

.47932 

-65.21 

.99 

.34627 

-176.89 

.45 

.47809 

-66.76 

1.00 

.33333 

0.00 

.46 

.47683 

-68.32 

1.01 

.31967 

177.01 

.47 

.47553 

-69.88 

1.02 

.30552 

174.18 

.48 

.47419 

-71 .44 

1.03 

.29111 

171.53 

.49 

.47281 

-73.01 

1.04 

.27669 

169.08 

.50 

.47140 

-74.58 

1.05 

.26247 

166.85 

.51 

.46996 

-76.15 

1.06 

.24863 

164.85 

.52 

.46849 

-77.73 

1.07 

.23534 

163.08 

.53 

.46699 

-79.31 

1.08 

.22270 

161.54 

.54 

.46547 

-80.89 

1.09 

.21078 

160.21 

.55 

.4639? 

-82.48 

1.10 

.19963 

159.08 

.56 

.46236 

-84.07 

1.11 

.18927 

158.15 

.57 

.46077 

-85.66 

1.12 

.17969 

157.39 

.58 

.45918 

-87.26 

1.13 

.17086 

156.78 

.59 

.45757 

-88.86 

1.14 

.16275 

156.32 

.60 

.45596 

-90.47 

1.15 

.15531 

155.97 

.61 

.45434 

-92.08 

1.16 

.14851 

155.72 

.62 

.45273 

-93.69 

1.17 

.14229 

155.55 

.63 

.45113 

-95.32 

1.18 

.13659 

155.46 

.64 

.44954 

-96.95 

1.19 

.13138 

155.41 
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C   C   VTF(T(5)-PHASF)  VS.  FREQUENCY.  OUTPUT  DATA,  TYPE  2  FILTER,  FIG. 31 


Wl 


T(  JW) 


Wl 


T(JW) 


01 


.10 

.50000 

-15.01 

.65 

.49079 

-107.20 

.]  1 

.50000 

-16.52 

.66 

.48988 

-109.24 

.12 

.50000 

-18.03 

.67 

.48888 

-111.31 

.14 

.50000 

-21.06 

.69 

.48660 

-115.51 

.15 

.49999 

-22.57 

.70 

.48529 

-117.66 

.17 

.49999 

-25.61 

.72 

.48230 

-122.04 

.19 

.49998 

-28.66 

.74 

.47871 

-126.54 

.  2  0 

.49998 

-30.19 

.75 

.47666 

-128.84 

.22 

.49997 

-33.27 

.77 

.47194 

-133.55 

.23 

.49996 

-34.81 

.78 

.46923 

-135.96 

.24 

.49995 

-36.35 

.79 

.46627 

-138.41 

.25 

.49994 

-37.90 

.80 

.46303 

-140.90 

.26 

.49993 

-39.46 

.81 

.45948 

-143.42 

.27 

.49992 

-41.02 

.82 

.45561 

-145.98 

.28 

•49990 

-42.58 

.83 

.45139 

-148.58 

.29 

.49989 

-44.15 

.84 

.44678 

-151.21 

.30 

.49986 

-45.73 

.85 

.44178 

-153.88 

.31 

.49984 

-47.31 

.86 

.43636 

-156.59 

.32 

.49981 

-48.90 

.87 

.43049 

-159.32 

.33 

.49978 

-50.49 

.88 

.42416 

-162.08 

.34 

.49975 

-52.09 

.89 

.41736 

-164.87 

.35 

.49970 

-53.70 

.90 

.41009 

-167.68 

.36 

.49966 

-55.31 

.91 

.40233 

-170.50 

.37 

.49961 

-56.93 

.92 

.39410 

-173.33 

.38 

.49955 

-58.56 

.93 

.38540 

-176.16 

.39 

.49948 

-60.20 

.94 

.37627 

-178.99 

.40 

.49941 

-61.85 

.95 

.36672 

-181.81 

.41 

.49932 

-63.50 

.96 

.35679 

-184.61 

.42 

.49923 

-65.32 

.97 

.34653 

-187.38 

.43 

.49913 

-67.00 

.98 

.33599 

-190.11 

.44 

.49901 

-68.68 

.99 

.32522 

-192.80 

.45 

.49888 

-70.37 

1.00 

.31427 

74.48 

.46 

.49874 

-72.08 

1.01 

.30322 

161.83 

.47 

.49858 

-73.80 

1.02 

.29211 

159.32 

.48 

.49841 

-75.52 

1.03 

.28102 

156.88 

.49 

.49822 

-77.26 

1.04 

.26999 

154.53 

.50 

.49800 

-79.01 

1.05 

.25909 

152.25 

.51 

.49777 

-80.78 

1.06 

.24837 

150.07 

.52 

.49751 

-82.56 

1.07 

.23786 

147.98 

.53 

.49723 

-84.35 

1.08 

.22760 

145.98 

.54 

.49692 

-86.16 

1.09 

.21764 

144.08 

.55 

.49658 

-37.98 

1.10 

.20799 

142.27 

.56 

.496  20 

-39.82 

1.11 

.19868 

140.57 

.57 

.49579 

-91.68 

1.12 

.18973 

139.12 

.58 

.49535 

-93.55 

1.13 

.18114 

137.60 

.59 

.49486 

-95.44 

1.14 

.17291 

136.19 

.60 

.49432 

-97.35 

1.15 

.16507 

134.87 

.61 

.49373 

-99.28 

1.16 

.15759 

133.64 

.62 

.49309 

-101.23 

1.17 

.15047 

132.50 

.63 

.49239 

-103.20 

1.18 

.14372 

131.44 
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C   VTF (T(S)-PHASE)  VS.  FREQUENCY.  OUTPUT  DATA,  TYPE  3  FILTER.  FIG. 34 


Wl 


T(  JW) 


01 


Wl 


T  (  JW) 


01 


.10 

•  50000 

-15.67 

.65 

.48203 

-113.19 

.11 

.50000 

-17.25 

.66 

.48023 

-115.37 

.12 

.50000 

-18.82 

.67 

.47828 

-117.58 

.14 

•  50000 

-21.98 

.69 

.47390 

-122.06 

.15 

.50000 

-23.57 

.70 

.47144 

-124.34 

.17 

.49999 

-26.74 

.72 

.46594 

-128.98 

.19 

.49999 

-29.93 

.74 

.45957 

-133.72 

.20 

.49999 

-31.53 

.75 

.45604 

-136.13 

.22 

.49998 

-34.75 

.77 

.44821 

-141.02 

.23 

.49997 

-36.36 

.78 

.44389 

-143.50 

.24 

.49996 

-37.98 

.79 

.43929 

-146.01 

.25 

.49995 

-39.61 

.80 

.43440 

-148.54 

.26 

.49993 

-41.23 

.81 

.42921 

-151.09 

.27 

.49992 

-42.87 

.82 

.42371 

-153.66 

.28 

•49990 

-44.5] 

.83 

.41789 

-156.25 

.29 

.49987 

-46.16 

.84 

.41176 

-158.85 

.30 

.49984 

-47.81 

.85 

.40531 

-161 .47 

.31 

.49981 

-49.47 

.86 

.39854 

-164.10 

.32 

.49977 

-51.14 

.87 

.39145 

-166.74 

.33 

.49972 

-52.82 

.88 

.38406 

-169.38 

.34 

.49967 

-54.50 

.89 

.37638 

-172.03 

.35 

.49960 

-56.19 

.90 

.36841 

-174.67 

.36 

.49953 

-57.90 

.91 

.36017 

-177.31 

.37 

.49944 

-59.61 

.92 

.35168 

-179.94 

.38 

.49934 

-61.33 

.93 

.34296 

-182.56 

.39 

.49923 

-63.05 

.94 

.33404 

-185.15 

.40 

.49910 

-64.79 

.95 

.32495 

-187.73 

.41 

.49895 

-66.70 

.96 

.31572 

-190.27 

.42 

.49879 

-68.46 

.97 

.30638 

-192.79 

.43 

.49860 

-70.24 

.98 

.29695 

-195.26 

.44 

.49839 

-72.02 

.99 

.28748 

-197.70 

.45 

.49815 

-73.82 

1.00 

.27800 

-200.09 

.46 

.49788 

-75.63 

1.01 

.26854 

157.41 

.47 

.49756 

-77.46 

1.02 

.25913 

155.12 

.48 

.49725 

-79.29 

1.03 

.24980 

152.89 

.49 

.49688 

-81.15 

1.04 

.24059 

150.71 

.50 

.49646 

-83.01 

1.05 

.23152 

148.60 

.51 

.49601 

-34. 9C 

1.06 

.22262 

146.55 

.52 

.49550 

-86.80 

1.07 

.21390 

144,57 

.53 

.49494 

-88.71 

1.08 

.20539 

142.65 

.54 

.49432 

-90.64 

1.09 

.19710 

140.81 

.55 

.49363 

-92.60 

1.10 

.18904 

139.04 

.56 

.49288 

-94.56 

1.11 

.18124 

137.34 

.57 

.49206 

-96.55 

1.12 

.17369 

135.71 

.58 

.49115 

-98.56 

1.13 

.16640 

134.15 

.59 

.49016 

-100.58 

1.14 

.15938 

132.66 

.60 

.48908 

-102.63 

1.15 

.15263 
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C   C   VTF(T(S)-PHASE)  VS.  FREQUENCY*  CUTPUl  DATA,  TYPE  4  FILTER,  FIG. 37 
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C   VTF(T(S)-PHASE)  VS.  FREQUENCY,  OUTPUT  DATA,  TYPE  6  FILTER,  FIG. 42 
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C   C   VTF(KS)-PHASF)  VS.  FREQUENCY*  OUTPUT  DATA,  TYPE  7  FILTER,  FIG. 45 
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APPENDIX  I 

Complementary  Filter  ADPI  Coefficient  Symmetry.   It 
can  be  noted  that  complementary  filter  ADPI ' s  exhibit  coef- 
ficient symmetry  in  their  numerator  and  denominator  poly- 
nomials.  This  observation  can  be  verified  mathematically. 

A  polynomial  P(s)  of  degree  n  has  coefficient  symmetry 
if 

snP(l/s)  =  P(s)  (50) 

Let  Z(s)=P(s)/Q(s)  be  the  driving  point  impedance  of  the 
low-pass  filter  of  the  complementary  filter  and  let  deg  P=m, 
deg  Q=m-1.   Then  the  high-pass  filter  has  a  driving  point 
impedance  Z(l/s)  because  the  low-pass  to  high-pass  trans- 
formation is  s  replaced  by  1/s. 

The  driving  point  impedance  of  the  complementary  filter 


is 


«(.)  „   Z(s)  Z<Vs)  (51) 


Z(s)  +  Z(l/s) 


Z  (s)  =   r(s)  P(Vs) (52) 

P(s)  Q(l/s)  +  P(l/s)  Q(s) 


z   (S)  -   P(s)   sroP(l/s) (53) 

CF       sm  P(s)  QCl/s)   +  sm  PCl/s)  QCs) 


The  test  for  coefficient  symmetry  can  now  be  applied  to 
the  numerator  first  and  then  to  the  denominator  of  Zcp(s). 
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The  numerator  degree  is  2m  and  since 


s2m[p(S)smP(l/s)] 


m 


s->l/s 


=  su,P(l/s)P(s) 


(54) 


the  numerator  polynomial  is  symmetric.   The  degree  of  the 
denominator  is  eaual  to  2m.   Hence 

s2m[smP(s)0(l/s)  +  smP(l/s)Q(s)l 

s-*l/s 

s2m[s-mP(l/s)Q(s)  ♦  s~mP(s)Q(l/s)]  =  (55) 

smP(s)Q(l/s)  +  smP(l/s;Q(s) 


implies  that  the  denominator  polynomial  is  symmetric. 

Note  that  the  numerator  and  denominator  polynomials 
degrees  are  both  equal  to  2m. 
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APPENDIX  J 

FORTRAN  Program  for  Coefficient  Evaluation  of  the  Pro- 
duct of  Two  Polynomials  with  Algebraic  Coefficients.   This 
computer  program  is  a  modification  of  one  written  by  Dr. 
Benton  Weathers.   Since  this  program  is  quite  long  end  rather 
complicated,  a  self-explanatory  program  was  added  instead 
of  a  block  diagram.   The  idea  behind  the  program  is  to  set 
up  two  three-dimensional  matrices  with  each  matrix  repre- 
senting one  of  the  polynomials.   Since  each  term  of  the 

■z     z  p 

polynomial  has  the  form,  8a  b;s  ,  each  dimension  of  the 
matrix  represents  an  exponent  of  a,  b,  or  s.   And  the  num- 
erical coefficient  was  assigned  to  the  location  given  by 
the  exponents.   Having  established  this,  the  multiplying 
and  term  collecting  operations  are  next. 
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APPENDIX  K 

ADPI  Parameter  Evaluation  for  a  Type  £  Filter.   In 
applying  the  aidentity  algorithm  to  the  ADPI  of  Type  7, 
three  equations,  (2$),  (24),  and  (25),  and  three  unknowns, 
C,  C-.,  and  m,  result.   Equation  (23)  can  be  solved  for  C-. 
in  terms  of  C  and  m.   This  results  in  equation  (56)  wh-re 
a=l+m.   Equation  (24)  can  he  manipulated  to  give  equation 


01   =   P  p     aC (56) 

1  a^CT  +  aC  -  1 


(57)    and   if   the   expression  for  C-,    in  equation   (56)    is   sub- 
stituted  for  the  C-.    in   equation   (57),    the   expression   for 
C   given   in   equation    (58)   will   be   obtained. 

C1(a5bCm+a5bG2m+a?C?+a2C-a2bC-a2b-aGm-l) 

+   a5bC2   -   a5bC   -  aC    =   0  (57) 


C    =  ^3  -  m2/(l+m)2  (58) 


Equation  (25)  can  be  rearranged  to  give  equation  (59). 
Substituting  the  expression  for  C-,  from  equation  (56)  into 
equation  (59)  results  in  equation  (61)  after  much  expanding, 
manipulating,  consolidating,  and  factoring.   Now  since  equa- 
tion (56)  has  been  used  with  both  (57)  and  (59)  to  give 
(58)  and  (61),  these  two  will  now  be  combined  to  give  a 

polynomial  with  only  one  parameter,  m.   Using  equatiqn  (58) 

4    5       2 

C  ,  C  ,  and  C   can  be  generated.   Substituting  these  into 
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C12X   +   C^Y  +    aCZ    =   0  '  (59) 


where 


X   =  C2a2b  +   2Ca2b  -  Cabm  +   a  b  -  abm  -  b 

Y  =  C2a4b2  +  C2a5  +  Ca4b2  +  Ca2m2  -  ab  -  a2b  -  ab        (60) 

-   a  -  1 
Z  =  Ca^b  +  Ca  -   a  b  -   ab  -   a  -  1 


AC4  +   BC5  +  DC2  +   EC   +    F   =   0  (61) 


where  A   =    (2+m-m  -nr  )(l+m) 

B  =    (2+m-2m2-m5+m4)(l+m)5 

D   =    (-7-3m+2m2+2m5+2m4+m5)(l+m)2  (62) 

E  b    (-4-2m+2m2+m5+m^)(l+ra) 

p^45 
F   =   6+2m-3m  -m   -m  -nr 

equation  (61)  and  clearing  denominators  results  in  equation 

(63)  after  factoring  and  consolidating.   Equation  (63)  can 

not  be  factored  for  roots  of  m.   Therefore  notice  that  eauation 

(3+2m-llm2-7ra5+llm4+7m5-3m6-2m7)  +  (2+m-6m2 
-3m5+5m4+2m5-m6)  I J   -  m2  =  0  (63) 

(63)  has  the  form  J>+QsJr=0.      If  this  is  multiplied  by  the 
factor,  (P-Q^/r),  the  result  will  be  P  -Q  R=0  which  can  be 
factored;  however,  this  polynomial  now  has  in  it  spurious 
roots  of  the  factor,  (P-C^/R) .   So  any  roots  obtained  will 
have  to  be  checked  in  equation  (  S3   )  to  make  sure  that  they 
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are  not  just  roots  of  the  factor,  CP-QfJR) .      Hofman's  computer 
program  of  Appendix  M  and  N  using  BAIRSTOW'S  method  on  the 
IBM  14-01-1410  gave  the  roots  to  the  fourteenth  order  poly- 
nomial.  These  roots  are  given  in  equation  (64)  and  only 
one  was  realizable  and  not  trivial.   Root  q,.  is  the  only 
realizable  non-trivial  root  and  when  substituted  into  equation 
(63)  is  found  to  be  a  root.   With  this  value  of  m  substituted 
into  equation  (58),  C  is  found  to  be  a   realizable  number  when 

14 
P2  -  Q2R  =  / |  (m+qn)  (64) 

n=l 

where  q,  =   1  qg  =  0.88  -  jO.49 

q2  =   1  Qo  =  0-88  +  J0.4Q 

Q3  =  -1  Q10  =  -1.59 

q4  =  -1  qn  =  -0.245  -  jO.515     (65) 

q^  =   0.7264  q12  =  -0.245  +  JO. 515 

q6  =  1.7  q15  =  -1.594 

q?  =   1.25  q1Z|.  =  -0.7143309 

evaluated  and  then  C-,  is  evaluated  from  equation  (56). 
This  gives 

m  =  0.7143309 

C   =  0.92041  (66) 

C1  =  0.51437 
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APPENDIX  L 

FORTRAN  Program  for  Evaluation  of  Type  7  Filter  ADPI 
Coefficients.   With  th«*  parameter  values  of  equation  (66), 
the  coefficients  expressed  algebraically  in  eauations  (22) 
were  evaluated  with  the  computer  program  given  on  pages 
126  and  127  with  the  results  given  on  page  127. 
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»F10«5»5X» 10H  B 
=  ,F10.5,5X,3H 


MON$$      JOB   FILTER  COEFFICIENTS  CASE 

MCN$$      CCMT  15  MIN,  10  PAGES*  FOWLER 

MON$$      ASGM  MJR»12 

MON$$      ASGN  MG0»16 

MON$$      MODE  GO»TEST 

MCN$$      EXEQ  FORTRAN, ,,,,, ,FANOUT 

DIMENSION  C(  10) ,D( 10 ) 
9  E0RMAT(5X,10H  A(  0) 
2  FORMAT(5X,3H  A( ,I2»5H J 

FM=0. 7143309 

F<1=0. 92041056 

FK2=0. 51436894 

A=1.+FM 

B=1.-FM 

AB=A*B 

A2B=AB*4 

A3B=A2B*A 

A4B=A3B*A 

A5B=A4B*A 

A6B2=A5B*A*B 

A2=A*A 

A3=A2*A 

A4=A3*A 

BK2=B*F<2 

FK1K?=FK1*FK2 

FK22=FK2*FK2 

FKll=FKl*FKl 

F1=1.+FK1 

F2=l.+FKl*FM 

F3=1.+BK2 

F4=1.+FK2*FM 

C0=A3B*FK1K2 

C( 1  )=A4B*FK1K2*F1+A*A*FK1K2 

C(2  )=A4B*FK1*(FK1+FK2*F2)+   A3*FK1K2*F 

C1=A5B*FK1*(FK1+BK2*F1 ) +A3*FK 1* ( FK 1+FK 

C(3 )=C1+A*FK2*(FK1+FK2*F2) 

C?  =  A5B*FKn+A4*FKl*(FKl+BK2*Fl)+A3*(FK 

C(4)^C2+A2*FK2*F1*(FK1+FK2*F2)+A*FK22* 

C3=A6B2*FK11+A4*FK11+A4*(FK1+BK2*F1 )** 

C(5)=C3+A2*FK22*F1**2 

C(6)=C(4) 

C(7)=C(3) 

C(8  )=C(2) 

C(9)=C(1) 

C(1C)=C0 

DO  =  CO 

D( 1  )=A3B*FK1*F4  +  A2*FK1K2  +  A2*BK2 

D(2  )=A4B*FK1*F3+A2*FK1*F4+A2*FK2*( FK1+ 

Dl=A4B*F(a+A3*FKl*F3+A2*F4*<FKl+BK2*Fl 

D(3  )=D1+A3B*(FK1+FK2*F2)+A2*FK2*F1+A*F 

D2=A5B*B*FK1+A3*FK1+A3*F3*(FK1+BK2*F1) 

D3=A*FK2  2*F1+A4B*(FK1+BK2*F1 )+A2*( FK1+ 

D(4)=D2+D3+FK2*F4 

D4  =  A4B*,rKl+A3*(  FK1+BK2*F1  )+A2*F3*(  FK1  + 


7 
♦  EE 


23  DEC  64 


(  C)   =  ,F10.5) 
B(»I2,5H)   =  ,F10.5) 


1+   A2*FK2*(FK1+BK2*F1 ) 
2*F2)+A3*FK2*Fl*( FKl+BK2*Fl ) 

1+BK2*F] )*(FK1+FK2*F2) 

Fl 

2+A2*(FKl+FK2*F2)**2+FK22 


BK2*F1 )+A3B*FK2*Fl+A*FK2 
)+A*FK2*( FK1+FK2*F2 ) 
K2*F3 

+A*F4*(FK1+FK2*F2) 
FK2*F2 ) +A2*FK2*F1*F3 

FK2*F2)+A*FK2*F1*F4+FK22 
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D(5  )=2.*D4 

D(6)=D(4) 

D(7 )=D(3) 

D(8 )=D(2) 

D(9)=D(1 ) 

D( 10)=D0 

WRITE  (2t9)C0,D0 

DC  3  1=1,10 

WRITE! 2»2) I »C(  I  )» I  »D(  I  > 
1  CONTINUE 

STOP 

END 
MCN$$      EXEQ  LINKLCAD 

CALL  FANCUT 
MCNSJ      EXEQ  FANCUT.MJB 


KCN$$  SL'y    JOB   FILTER  COEFFICIENTS  CASE  7 


23  DEC  64 


A( 

0) 

A  ( 

1  ) 

A( 

2) 

A( 

3) 

A( 

4) 

A( 

5) 

A( 

6) 

M 

7) 

A  ( 

8) 

A( 

9) 

A(10) 


.68140 

3,63470 

10.42514 

20.45234 

29.90444 

33.9341  1 

29.90444 

20.45234 

10.42514 

3.63470 

.68140 


B( 

0) 

B( 

1  ) 

B( 

2  ) 

B( 

3) 

B( 

4) 

B( 

5) 

B< 

6) 

B( 

7) 

B( 

8  ) 

B( 

9) 

B(  10) 


.68140 

3.63470 

10.42514 

20.45234 

29.93704 

33.77306 

29.93704 

20.45234 

10*42514 

3.63470 

.68140 
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APPENDIX  M 

Bairstow' 3  Root  Extraction.   Mr.  Larry  Hofman  prepared 
the  following  analysis  of  root  extraction  using  Bairstow1 s 
method  and  Hamming's  (12)  suggested  numerical  methods. 
However  Hofman  added  procedures  to  insure  convergence  for 
most  polynomials. 

Bairstow' s  Method.   Let  the  polynomial  to  be  factored  be 


P(x)  =  anxn  +  an_i: 


Ji-l 


+  a2X  +  a-^x  +  aQ  (67) 


There  is  a  quadratic  factor  of  the  form  Ax  +  Bx  +  C.   Assume 

a  =  1  and  guess  at  the  factor 
n         " 


x  +  px  +  q 


(68) 


Divide  the  polynomial  by  the  quadratic  factor  and  obtain 
a  quotient  and  a  remainder,  e.g., 


ax  +  a„ 
n      n-1 


n  1  P 

+  ...  +  a,x  +  an  =  (x  +  px 


l0 


n-2 


n-3 


+  qXb^11"*  +  bn_1xn  2  +  ...  +  b2)  +  bxx  +  bQ   (69) 


The  peculiar  subscripts  on  the  b's  make  notation  easier. 
In  a  skeleton  synthetic  form 


1  p  q 


an 

an-l 

an-2 

<3bn 

Pbn-1 

an-3 
"Vl 

Pbn-2 

... 
... 
... 

a2 

p^5 

al 
pb2 

a0 
qb2 



Pbn 

bn 

Vl 

bn-2 

v5 

... 

b2 

bl 

bo 

l?9 


where  the  remainder  is 

b,x  +  bQ 
The  algebraic  relations  between  the  coefficients  are 

bn   =  an 

Vl  =  an-l  "  ?bn 

bn-2  =  an-2  "  pVl  "  ^bn 

(70) 

bn-k  "  an-k  "  Pbn-k+l  "  aPn-k+2(k=2'  5'  ••••  n-1} 


bQ   =  s0  -  qb2 

The  desired  auadratic  factor  is  obtained  if,  and  only 
if,  the  remainder  is  identically  zero.   That  is, 

bl  -  bo  ■  ° 

Now  find  some  exact  method  of  correcting  the  guess  of  B 
and  C  in  order  that  the  above  conditions  are  met.   Consider 
the  coefficients  b,  and  bQ  to  be  functions  of  p  and  q: 

b0  =  bo^p,CI^ 

Using  Newton's  method  in  two  dimensions,  expand  b,  and  b^ 
about  the  present  guess  (p  and  a).   Writing  B  and  C  as  the 
desired  solution, 
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d  b,     a  b, 
^(B.C)  -  0  =  b1(p,q)+-g-^AP+-^7iAa+  •  ••     (71) 


8  b     e>bn 
bQ(B,C)  =  0  =  b0(p,o)+-^AP^AQ+ 


where  AP  =  B  -  p  (72) 

Aq  =  C  -  q 

are  the  errors  to  be  corrected  (approximately)  for  the  next 
guess.   Neglecting  all  but  the  linear  terms  in  (71),  results 
in  a  pair  of  linear  eouations  for  the  changes  to  be  made 
in  p  and  a.   The  problem  is  to  find  the  partial  derivatives 
which  are  the  coefficients  of  the  unknowns  AP  and  A°» 
Differentiating  eouations  (70)  with  respect  to  p,  gives 


3b»      .0 


8P 


3V 

dp 

•1 

= 

~bn 

-  P 

©P 

<9P 

-2 

= 

"bn- 

-1 

-  P 

c)P 

-i 

c>bn 

dbn- 

<9P 

-k 

B 

"bn- 

-k+1 

-   P 

0\. 

ap 

-k+1 

-  p  ap 

-k+2 

ab0 

©p 

— 

ab2 

-qar 
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Now  write 


3*i 
ap 


■  -  c, 


(73) 


then 


n 


'n-1 


'n-2 


=  0 
=  b. 
=  b 


n 


n-1 


-pc 


n 
Pcn-1 


qc* 

A   n 


(74) 


cn-k  ~  bn-k+l  "  pcn-k+l  -  qcn-k+2 


'0 


-  qc 


2* 


These  cauations  are  practically  in  the  same  form  as  (70). 

This  suggests  repeating  the  process  of  synthetic  division 

2 
using  the  quadratic  factor  of  x  +  px  +  q  on  the  b's  to 

obtain  coefficients  c^.  Then 


1   p   q 


>     b  ,   b  0  b  , 
n     n-1   n-2   n-3 


b2   bx 


'0 


qcn   qcn-l  ' ' '   qc4-  qc3   qc2 


—    pcn   Pcn«i  pcn-2  *••   pc3  pc; 


where 


cn    cn-l  cn-2  cn-3 


c    *  b 
n      n 


'0* 


n-i    n-l   r   n 

c„  o  -  bv  0  -  pc  t    -  ac 

n-2    b-2   '  n-1      -  n 


(75) 
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cn-k  =  Vk  "  pcn-k+l  "  qcn-k+2 


(75) 


cl   =  bl   "  pc2     "  qc3 
cQ   =  bQ  -  Qc2 


The  partial  derivatives  desired  in  (71),  using  (73), 
are 

—  =  -c*  ^„  =  -c* 


3P  "  "1  8>p     "0 
Comparing  (74)  and  (75)  gives 

ck-l  =  ck     (k  =  n>  n"1'  •••  »  3»  2) 

c0   B  cl  "  bl  +  pc2 


Hence 


=  -c2  -~—   =  -(c1  -b1  +  pc2)   (76) 


ap 


Nov;  examine  the  process  for  the  partial  derivatives 
with  respect  to  q.   Again  differentiating  (70)  gives 


n   =  0 


£q 


=  0 


SV2 

£>q 

"bn 

3Vi 

pao 

©bn 

gVk     _b  QVk-fi       Q  Vk+2 

0o  "  n-k+2  ~paq  ~  q6>q 


d bo  ^  b2 

Nov;  set 


c)q      k 


to  get 


c»   .  0 


C**    s  0 

n-1    u 


cn-2  =  bn     -  pcnll   ~  qcn* 


n-k  •  n-k+2  -  pc**   ,   qcn-k+2 

r   n— k+1 
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6>\ 

=  -c**  (77) 


(78) 


c**   =  b2  -  qc|* 


Since  c**  =  c**..  =  0,  it  is  necessary  to  identify 

ck-2=  ck  ^k  =  n'  n"1'  '••  '  4'5^ 

c0*  =  c2  +  pc3 

if  (75)  and  (78)  are  to  be  compared. 

The  partial  derivatives  desired  for  (71)  are 
d)b.  £>b 


1    _..      '     v~0 


-C**  =  -Cz    -=c—  =  -c**  =  -(c0  +  PC,). 


-^-  -  -^  -  -„3    ^- 0  -  -vv,2  -r  v-3 
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Thus      0-^,0,  =  c2^P  +  c^  AQ  (79) 

t>0(p,q)  =  (c1  ~  b^  +   pc2)AP  +(c2+pc5)  An 

Solving  equation  (79)  simultaneously,  results  in 

b, (cp  +  pc.)  -  c^bn 

Ap  =  J^J. — p L2  (80) 

c2  +  c5(b1  -  C]_) 

b0c2  "  ^l^cl  "  bl  +  pc2^ 
AQ  =  — p 

c2  +  c3^bl  "  cl^ 

Now  replace  the  value  of  p  and  a  with 

p — o  P  +  AP 
a  — >  a  +  AQ 

and  repeat  the  above  process  until  the  values  of  b,  and 
bQ  are  sufficiently  small.   The  convergence,  when  it  works, 
is  quadratic;  that  is,  the  errors,  when  small,  are  approx- 
imately souared  each  step.   Thus  an  iterative  process  is 
obtained  which  will  converge  upon  a  Quadratic  factor  of  the 
original  polynomial.   V/hen  found,  it  can  be  factored  out 
(the  first  division  step)  and  then  use  the  quotient  as  a 
new  polynomial  to  be  examined  by  the  same  process. 

The  errors  introduced  by  this  method  are  accumulated 
as  each  factor  is  removed  (division  by  an  inexact  root  and 
then  discarding  the  remainder).   Thus  a  polynomial  of  large 
order  will  reauire  a  high  degree  of  accuracy  in  order  that 
the  last  factors  will  have  an  acceptable  accuracy  also. 
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The  convergence  of  Bairstow's  method  is  somewhat  sen- 
sitive to  the  initial  guess  of  p  and  o.   Since  no  value 
has  been  found  for  p  and  q  which  will  cause  convergence  for 
all  cases,  Hofman  devised  a  program  that  tries  a  succession 
of  initial  guesses,  as  outlined  in  Table  I,  in  an  attempt 
to  extract  a  stubborn  factor.   This  program  has  yet  to  fail 
after  many  trials. 


trial 

P 

a 

roots 

1 

4 

3 

-1   -3 

2 

2 

1 

-1   -1 

3 

0 

-1 

-1   +1 

4 

2 

2 

-l  ±  di 

5 

-2 

2 

+i  t  di 

Table  I 
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APPENDIX  N 

Hof man ' s  FORTRAN  Program  for  Polynomial  Root  Extraction, 
With  the  analysis  for  root  extraction  as  presented  in  Appendix 
M,  Hofman  wrote  the  following  FORTRAN  computer  program  for 
use  on  the  IBM  1401-1410  digital  computer. 
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PRINT 

2 

x  +px+o 


READ  POTT.  ORDER,  N 


READ  COEFFICIENTS,  ap 


SET  INITIAL  ESTIMATES 
FOR  p  and  a 


compute  c ' s 

s 

! 

compute 
p   and        q 

i 

P  =  P  +    AP 
a   =   a+    A° 

PRINT 
2 

X+PX+Q 


reset  a's 
from  b' s 


N  =  N-2 


Fig.  59.   Block  diagram  of  Bairstow's  polynomial  factoring 
computer  program. 
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1 
2 
3 
4 
5 
6 
7 
8 
9 

In 

C 


C 
1? 


c 

14 


C 
15 

16 

17 
18 

C 
20 


CF  POLYNOMIAL 
HOFMAN,  ELECT 


FACTOR 
ENGG 


1/1/64 


MON$$      JOB   BAIRSTOW  METHOD 

MONSS      COMT  5  MIN,  5  PAGES* 

MON$$      ASGN  MJB»12 

MON$$      ASGN  MG0»\6 

MON$$      MODE  GO 

MON$$      EXEQ  FORTRAN, , ,18,3, , ,BAIRSTOW 

DIMENSION  A(2C),B(20),C(20),PI(5),QI(5) 

EQUIVALENCE  (Al,A(l)),(Bl,B(l)),(B2,B(2)),(B3,B(3)),(Cl,C(l)), 
1(C2,C(2)),(C3,C(3)),(C4,C(4)) 

FORMAT( 12) 

F0RMAT(E14.8) 

FORMAT( 18HKMETH0D  HAS  FAILED) 

FORMAT( 1X,F16.8) 

F0RMAT(4HSX  =,F11.7,25H  » 


LBHJ 


+  J  >!■ 


1 1.7,3X»7H 


F0RMAT(4HSX  =,F11.7,4H 
F0RMATI4HSX  =,F11.7) 

F0RMAT(39HlC0EFFICIENTS  IN  DECREASING  POWERS 
FORMAT( 1HK,60X»61HITERATI0N     B(l) 
1  Q/) 

F0RMAT(61X,I5,3X,1PE14.7,3(1PE15.7)) 
SET  5  INITIAL  GUESSES, 
PI  (  1  )=4. 
QI  (  1  )=3. 
PI (2)=2. 
01 (2)=1. 
PI  (3)=0. 
QI <3)=-l. 
PI (4)=2. 
QI (4)=2. 
PI ( 5)=-2. 
QI ( 5)=2. 

READ  ORCIR  N.   N=0  TO  TERMINATE  PROGRAM, 
RFAD( 1 ,1) N 
IF(N.EQ.O)  STOP 
N  =  N  +  1 

WPITE(3*8) 
DO  14  I  =  1,N 
J=N+1-I 

READ  COEFFICIENTS  IN  DECREASING  POWERS. 
READ( 1,2) A( J) 
WRITE(3,4)A( J) 
IF( A(N) .EQ.l. )  GO  TO  16 
DO  15  1=1. N 

NORMALIZE  COEFFICIENTS. 
A(  I  )=A(  I ) /A(N) 
1=1 

WRITE(3,9) 
IFIN-3)  130,120,18 
L  =  N-2 
ITRY=1 

SET  INITIAL  GUESS  FOR  P 
P=PI ( ITRY) 
Q  =  QI  (  ITRY) 
ITCNT=1 
CALCULATE  BS. 


X  =,F11.7) 
=,F11.7,4H  - 

OF  X/  ) 
B(0) 


J,F11.7) 


AND  Q. 
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2R 


30 


B(N)  = 
B(N-1 


1 


4* 
C 


50 


5* 

C 
57 


C 

6o 


80 
90 

no 


C 
1?0 


c 

130 


)=A(N 
DC  30  K=2. 
M  =  N-K 

B(M)=A(M)- 
B1=A1-Q*B3 
WRITE(3,10 
CHECK  ACCU 
IF(ABS(B2) 
IF(ABS(B1) 
ITCNT=ITCN 
IF(  ITCNT.G 
CALCULATE 
C(N)=1. 
C(N-1)=B(N 
DC  50  K=2, 
M  =  N-K 
C(M)  = 
C1  =  B1 


-1)-P 
L 

P*B(M+1  )-Q*B(M-t-2  ) 

) ITCNT»B2.B1,P»Q 

RACY  CF  GUESS. 

•  GE. • 00000001 )  GC  TC  45 

. LT.. 00000001 )  GC  TC  60 

T+l 

Ti25)  GC  TC  150 

CS  FCR  CORRECTION  CF  P  AND  0. 

-D-P 
L 


B(M)-P*C(M+1 )-Q*C(M+2 ) 
-Q*C3 

DENCM=C3*C 


3+C4*(B2-C2 ) 
NCM.EQ.O. )  GC  TC  55 

DELTA  P  AND  DELTA  Q. 
(C3+P*C4)-C4*B1 ) /DFNOM 
B1-(B2*(C2-B2+P*C3) ) )/DENCM 


LAli 
=  (B2* 
=  (C3*I 
57 


IF(DF 

CALCU 

DFLTP 

DFLTO 

GC  TC 

DFLTP=.l 

DELTQ=.l 

CORRECT  P  AND  Q. 

P=P+DELTP 

Q=0+DELTQ 

GC  TC  25 

ROUTINE  T! 

DSCRM=P*P- 

IF(DSCRM.L 

RCCT1=(-P+ 

RCCT2=(-P- 

WRITF(3»5 )R0CT1 ,R0CT2 

N  =  N-2 

DO  90  J=ltN 

A ( J)=B( J+2) 

0  (  17,12)  > I 


FACTOR  QUADRATIC. 
4.*Q 

T.O. )  GO  TC  110 
SQRT(DSCRM) )*.5 
SQRT(DSCRM) )*.5 


GC  T: 
REAL 
CXPT 
WRIT 
GO  T 
FACT 
1=2 
P  =  A( 
0  =  A( 


=-P*.5 

=SQRT( 
E(3f6  ) 

0  80 

OR  LAST  QUADRATIC  REMAINING. 


-DSCRM)*.5 
REAL»CXPT,REAL»CXPT 


2) 

1  ) 


GC  TC  60 

REMC 

ROOT 


WRIT 
GC  T 


VE  LAST  LINEAR  FACTOR 
=  -A( 1  ) 
F(3»7)RCCT 
C  12 
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C      ppFPARF  FOR  ANOTHER  INITIAL  GUFSS 
150    ITRY=ITRY+1 

IF(  ITRY.LE.5 )  GO  TO  20 

WRITE(3»3 ) 

GO  TO  1? 


END 
MON$$ 

MON$$ 


EXEQ  LINKLOAD 
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Zobel,  Bode,  Guillemin,  Norton,  Rowlands,  and  Szentirmai 
used  several  methods  with  limited  success  to  improve  the 
driving  point  impedance  of  fan-out  filter  configurations. 
Fritzemeyer  used  these  ideas  and  King's  approximate  identity 
in  presenting  the  aidentity  algorithm  and  impedance  elision 
design  procedures  for  complementary  filters. 

This  thesis  presents  in  detail  these  two  design  proce- 
dures and  expands  on  the  analysis  of  Fritzemeyer' s  interactance 
(a  fan-out  filter  figure  of  merit).   The  Zobel  process  for 
improvement  of  characteristic  impedance  and  the  driving 
point  impedance  is  the  focus  of  interest. 

This  investigation  defines,  codifies,  and  verifies  the 
Zobel  process,  aidentity  algorithm,  and  impedance  elision 
8S  methods  of  improving  the  frequency  response  characteris- 
tics of  constant -k  complementary  filters.   It  is  shown  that 
these  methods  not  only  improve  the  aidentity  driving  point 
impedance  or  interactance,  but  also  the  voltage  transfer 
function  with  linear  phase  shift. 

Verification  is  shown  with  graphs  of  aidentity  driving 
point  impedance,  interactance,  and  voltage  transfer  function 
versus  frequency.   These  characteristics  were  computed  on 
the  IBM  1620  and  IBM  1401-1410  digital  computers  using  FORGO 
and  FORTRAN  languages. 

In  the  course  of  this  verification  several  related 
topics  of  interest  such  as  the  FORGO  program  for  evaluating 
the  coefficients  of  a  specialized  "ladder"  network's  VTF, 


the  FORTRAN  program  for  coefficient  evaluation  of  the  product 
of  two  polynomials  with  algebraic  coefficients,  the  proof 
of  the  complementary  filter  ADPI  coefficient  symmetry,  end 
the  ADPI  parameter  evaluation  for  a  Type  7  filter  were  studied 
and  are  included  in  the  appendices. 


